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ABSTRACT    #191 

Arsenic  contamination  of  surface  and  ground  waters  is  an  important  problem  in  many 
watersheds  in  Montana,  either  as  a  result  of  natural  processes  or  mining  activities.  The 
Madison  River  in  Southwestern  Montana  contains  elevated  As  concentrations  (100  to  300 
ug  L-1  )  as  a  result  of  natural  geothermal  sources  near  its  headwaters.  Previous  studies 
have  concluded  that  irrigation  with  Madison  River  water  may  be  responsible  for 
contamination  of  domestic  wells  near  Three  Forks,  MT.  The  objectives  of  this  study  were 
to  (I)  determine  sorption  capacities  of  irrigated  and  nonirrigated  soils  of  the  Madison  River 
basin,  (ii)  evaluate  the  potential  role  of  redox  reactions  on  the  solubility  and  subsequent 
transport  of  As  through  soils,  and  (iii)  determine  effects  of  pH  and  phosphate  on  As 
transport.  Sorption  experiments  conducted  under  batch  conditions  showed  that  Madison 
River  soils  with  long-term  irrigation  histories  had  sufficient  As  sorption  capacity  to  retain 
As  applied  in  irrigation  water.  Furthermore,  As  sorption  was  lowest  in  those  soils  which 
supported  a  lower  redox  status.  In  these  conditions,  a  significant  fraction  of  the  total 
soluble  As  was  present  as  arsenite  rather  than  arsenate.  Controlled  redox  experiments 
with  Madison  River  soils  and  representative  solid  phases  showed  that  the  solubility  of  As 
increases  as  redox  levels  approach  0  mV.  Column  experiments  showed  that  the  transport 
of  As  increased  with  increasing  pH  especially  at  pH  values  exceeding  8-8.5.  and  at  P:As 
molar  ratios  exceeding  10:1.  The  amount  of  As  transported  through  columns  containing  a 
representative  aquifer  sand  was  a  function  of  column  pore  water  velocity:  As  recovered  in 
the  effluent  increased  with  increasing  pore  water  velocity,  indicating  that  As  sorption 
reactions  were  kinetically  limited  during  transport.  Results  from  these  studies  indicate  that 
redox  conditions,  pH,  phosphate  concentration  and  pore  water  velocity  all  play  important 
roles  in  the  amount  of  As  transported  through  soils. 


December  1,  1999 

Copyright  Editor 

Journal  of  Environmental  Quality 
677  S.  Segoe  Road 
Madison,  Wl  53711 


■ 


M 


DEC  0  8  1999 


ONTANA 

University  System 


tf£ 


,    \%<& 


RE:       Permission  to  reprint/publish  3  articles 
Greetings! 


Water  Resources 

CENTER 


The  3  articles  listed  below  are  the  result  of  an  original  Montana  Water  Center  'seed'  research  project 
begun  in  FY94.   Based  on  the  success  of  this  initial  work,  the  P.l.s  were  successful  in  obtaining 
additional  funding  to  expand  their  project. 

Journal  of  Environmental  Quality 

Jones,  C.A.,,  Inskeep,  W.P.,  and  D.R.  Neuman,  1997,  Arsenic  Transport  in  Contaminated  Mine  Tailings 

following  Liming,  Vol.  26  no.  2,  March-April,  26:433-439. 

Darland,  J.E.,  and  W.P.  Inskeep,  1997,  Effects  of  pH  and  Phosphate  Competition  on  the  Transport  of 

Arsenate,  Vol.  26,  no..  4,  July-August,  26:1133-1139.  (acknowledges  Water  Center) 

Jones,  CA.,  Inskeep,  W.P.,  Bauder,  J.W.,  and  K.E.  Keith,  1999,  Arsenic  Solubility  and  Attenuation  in  Soils  of 

the  Madison  River  Basin,  Montana:  Impacts  of  Long-Term  Irrigation,  Vol.  28,  no.  4,  July  -  August,  28:1314- 

1320.  (acknowledges  Water  Center) 

An  initial  Water  Center  funding  request  is  that  P.l.s    "...  publish  their  findings  in  a  technical  journal  of 
their  choice."  As  part  of  its  mandate,  the  Center  also  includes  a  copy  of  these  findings  in  its  Technical 
Reports  Library'  -  a  small  collection  maintained  here  at  the  Center  and  at  Montana  university  libraries. 
What  I  need  from  you  is  your  permission  to  publish  and/or  reprint  these  articles  as  part  of 


Montana  University  System  Water  Center  Technical  Report  Series 
of  Arsenic  in  Soils  of  the  Madison  River  Basin. 


#191-  Fate  and  Mobility 


An  average  of  30  photocopies  will  be  made  for  distribution  free  of  charge.   Typical  users  of  our 
Technical  Report  Library  are  university  faculty/staff,  agency  professionals,  scientists,  and  engineers 
from  throughout  Montana.   FYI  -  abstracts  of  most  of  our  reports  are  available  on  the  WEB  at 
http://water.montana.edu/docs/reports.html 

Could  you  please  send  the  necessary  'permissions  to  publish'  forms?  Thanks  for  taking  the  time  on  this. 
Please  call  or  e-mail  me  if  you  have  any  further  questions. 

Sincerely,  r\ 


Kathy  Stephens,  Section  104(B)  Program  Manager 

ph:  406-994-1772 

Visit  MONTANA  WATER  at    http://water.montana.edu 

cc:  Dr.  William  Inskeep 


stephens@montana.edu 

C:\Main\Pubs\inskeep-jrnl-env.quality.wpd 


Dorothy  Bradley,  Director 

Huffman  Building  101 

P.O.  172690 

Montana  State  University 

Bozeman,  MT  59717-2690 

Fox  406-994-1774 

Tel.  406-994-6690 

E-mail   wwwrc@gemini.oscs.montana.edu 

Web  Site   http://www.montana.edu/wwwrc 


Marvin  Miller,  Assoc.  Director 

Montana  Bureau  of  Mines 

Montana  Tech 

1300  W.Park  Ave. 

Butte,  MT  59701-8997 

Fax  406-496-4451 

Tel.  406-496-4155 

E-mail   marv@mbmgsun.mtech.edu 


Don  Potts,  Assoc.  Director 

School  ot  Forestry 

University  of  Montana 

Missoula,  MT  59812 

Fax  406-243-4510 

Tel.  406-243-5522 

E-mail   fo_dfp@selway.umt.edu 


Reprinted  from  (lie  Journal  <>/  Environmental  Quality 

Volume  26,  no.  2,  March-April  1997,  Copyright  O  1997,  ASA.  CSSA.  SSSA 

677  South  Segoe  Rd  ,  Madison.  Wl  5371 1  USA 

Arsenic  Transport  in  Contaminated  Mine  Tailings  following  Liming 

C.  A.  Jones,  W.  P.  Inskeep,*  and  D.  R.  Neuman 


ABSTRACT 

The  practice  of  liming  to  remediate  contaminated  soils  and  mine 
tailings  has  the  potential  to  mobilize  arsenic  (As),  due  to  the  pH 
dependence  of  As  sorption  reactions  on  oxide  minerals  and  layer 
silicates.  The  objectives  of  this  study  were  to  determine  the  effects  of 
liming  on  As  mobility  in  mine  tailings  and  identify  possible  mechanisms 
controlling  As  mobilization  with  increased  pH.  Six  mine  tailing  samples 
obtained  from  an  abandoned  copper  smelter  near  Anaconda,  Montana 
were  analyzed  for  total  As  and  soluble  constituents  using  saturated 
paste  extractions.  Concentrations  of  soluble  As  among  the  six  samples 
did  not  correlate  with  total  As,  but  were  more  closely  related  to  pH. 
Saturation  indices  with  respect  to  known  metal  arsenate  solid  phases 
suggested  that  metal  arsenate  solid  phases  were  not  controlling  soluble 
As  in  these  samples.  Two  low  pH  samples  (pond  tailings  and  reproc- 
essed tailings)  were  chosen  for  more  detailed  chemical  characterization 
and  unsaturated  column  transport  experiments  before  and  after  lim- 
ing. Soluble  As  concentrations  measured  in  column  effluent  increased 
by  factors  of  10  (reprocessed  tailings,  RT)  to  400  (pond  tailings,  PT) 
following  liming.  Sequential  extractions  of  these  tailing  samples  showed 
that  the  PT  contained  significantly  higher  "labile"  As  relative  to  the 
RT,  consistent  with  amounts  of  As  mobilized  after  liming.  Further 
characterization  of  these  samples  using  scanning  electron  microscopy 
(SEM)  and  energy  dispersive  analysis  of  x-rays  (EDAX)  suggested  that 
these  samples  do  not  contain  discrete  metal  arsenate  solid  phases. 
Based  on  this  suite  of  experimental  data,  increased  mobility  of  As 
with  liming  appears  to  be  consistent  with  the  pH  dependence  of  sorption 
reactions  of  As  on  Fe  oxide  minerals  rather  than  dissolution- 
precipitation  reactions  involving  As. 


Liming  is  an  established,  effective  method  for  immobi- 
j  lizing  trace  metals,  elevating  pH,  and  p'  omoting 
plant  reestablishment  in  low  pH,  contaminated  mine 
tailings  (Neuman  et  al.,  1993b);  however,  increases  in 
soil  pH  may  result  in  As  mobilization  because  of  the 
pH  dependence  of  sorption  reactions.  This  raises  the 
concern  of  contaminating  drinking  water  supplies  follow- 
ing liming  in  soils  or  tailings  containing  elevated  levels 
of  As  that  are  hydraulically  connected  to  groundwater 
or  surface  water.  The  USEPA  may  reduce  the  maximum 
contaminant  level  for  As  in  drinking  water  from  0.67 
\iM  to  between  0.0067  and  0.134  uM,  due  to  recent 
evidence  that  the  current  criteria  may  result  in  human 
cancer  risks  as  high  as  13  in  1000  (Pontius  et  al.,  1994). 
Consequently,  there  is  significant  interest  in  processes 
controlling  the  transport  of  As  into  surface  or  ground- 
waters. 

Previous  studies  have  investigated  the  pH  effects  of 
As  sorption  on  soils  or  well-characterized  solid  phases 
(Hingston  et  al.,  1971;  Goldberg  and  Glaubig,  1988). 
Goldberg  and  Glaubig  (1988)  found  that  the  sorption  of 
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arsenate  by  montmorillonite  and  kaolinite  was  highest 
near  pH  5  and  dropped  significantly  between  pH  7  and  9. 
Hingston  et  al.  (1971)  determined  that  arsenate  sorption 
on  both  goethite  and  gibbsite  decreased  with  increasing 
pH.  Goldberg  (1986)  reported  good  correlations  between 
experimental  and  modeled  arsenate  sorption  on  Fe  and 
Al  oxides  using  a  constant  capacitance  model,  which 
assumed  a  ligand  exchange  mechanism.  Jacobs  et  al. 
(1970)  showed  that  the  elimination  of  amorphous  Fe  and 
Al  by  treatment  with  oxalate  substantially  reduced  or 
eliminated  arsenate  sorption  capacity  in  24  Wisconsin 
soils,  implying  soil  arsenate  was  primarily  sorbed  to 
oxides  of  Fe  and  Al.  In  contrast,  arsenate  solid  phases 
are  generally  believed  to  be  nonexistent  in  surface  soils, 
due  to  their  relatively  high  solubility,  except  at  high 
arsenate  levels  (Robins,  1981;  Sadiq  et  al.,  1983;  Dove 
and  Rimstidt,  1985;  Nishimura  et  al.,  1987).  For  exam- 
ple, Dove  and  Rimstidt  (1985)  found  that  at  an  Fe  activity 
of  1  mA/,  arsenate  activities  needed  to  be  at  least  0.01  M 
for  scorodite  (FeAs04  ■  2H20)  to  be  stable,  which  is  much 
higher  than  typically  found  in  soil  porewater.  Thus,  the 
effects  of  pH  on  arsenate  solubility  in  soils  are  due 
primarily  to  surface  complexation  reactions  rather  than 
precipitation-dissolution  reactions . 

One  of  the  recommended  liming  strategies  for  maxi- 
mum metal  immobilization  in  mine  tailings  of  the  Clark 
Fork  River  Basin  of  Montana  consists  of  raising  the  soil 
pH  to  near  10,  using  a  mixture  of  Ca(OH)2  and  CaCOi 
and  allowing  the  soil  to  equilibrate  for  approximately 
1  yr  prior  to  plant  reestablishment.  During  equilibration, 
the  pH  generally  declines  to  values  near  8  (Schafer  et 
al.,  1989).  Field  research  at  contaminated  mining  sites 
in  the  Clark  Fork  Basin  has  shown  that  this  liming 
strategy  results  in  metal  immobilization,  yet  may  lead 
to  As  mobilization  (Neuman  et  al.,  1993a). 

The  objectives  of  the  present  study  were  to:  (i)  deter- 
mine the  extent  of  As  mobilization  in  low  pH  mine  tailings 
after  liming  and  (ii)  evaluate  probable  mechanisms  for 
increased  As  transport  after  liming.  To  accomplish  these 
objectives,  several  representative  mine  tailings  and 
metal-contaminated  soils  were  characterized  and  two  of 
the  tailings  were  subjected  to  transport  experiments  with 
and  without  liming. 

MATERIALS  AND  METHODS 

Sample  Collection  and  Characterization 

Six  samples  (0  to  0.30-m  depth)  were  collected  within  2  km 
of  an  abandoned  copper  smelter  located  in  Anaconda,  MT, 
and  included  pond  tailings  (PT),  smelter  stack  deposits  (SSD), 
and  reprocessed  tailings  (RT).   All  sampling  locadons  were 


Abbreviations:  SEM,  scanning  electron  microscopy,  EDAX,  energy  dis- 
persive analysis  of  x-rays;  SSD.  smelter  stack  deposits;  RT.  reprocessed 
tailings;  DOC.  dissolved  organic  carbon;  ICP,  inductively  coupled  plasma 
spectrometry,  IC,  ion  chromatography;  CB,  citrate  bicarbonate,  CDB, 
citrate  dilhionite  bicarbonate,  PT.  pond  tailings 
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within  a  site  that  was  placed  on  the  National  Priority  List  in 
1983  due  to  metal  contamination  following  96  yr  of  copper 
smelting  (Neuman  et  al.,  1993b).  Samples  were  collected  in 
triplicate  at  two  locations.  Concentrations  of  total  As  and  total 
cations  (Ca,  Cu,  Fe,  Mn,  Pb,  and  Zn)  were  determined  using 
x-ray  fluorescence  (Jones,  1982)  and  particle  size  analysis  was 
performed  using  the  hydrometer  method  (Gee  and  Bauder, 
1986).  Saturated  pastes  (Rhoades,  1982)  were  centrifuged  at 
10  000  £  and  filtered  through  0.20-um  nylon  filters  to  character- 
ize soil  solution  chemical  content.  The  resulting  filtrates  were 
analyzed  for  dissolved  organic  carbon  (DOC)  using  a  Dohr- 
mann  DC  190;  cations,  B,  and  S  by  inductively  coupled  plasma 
spectrometry  (ICP);  and  CI  ,  NCK,  SOl  ,  and  POl~  by  ion 
chromatography  (IC).  Total  alkalinity  was  measured  by  HCI 
titration  to  an  inflection  point  near  pH  4.5.  These  titrated 
samples  were  then  sparged  with  N2  (g)  for  more  than  15  min, 
titrated  to  the  original  pH  with  NaOH,  and  retitrated  to  the 
inflection  pH  with  HCI  to  calculate  carbonate  and  noncarbonate 
alkalinity.  Continuous-flow  hydride  generation  atomic  absorp- 
tion spectrophotometry  (HG-AAS)  was  used  to  analyze  As. 
Samples  were  acidified  to  3M  HCI,  prereduced  with  1%  KI, 
and  treated  with  both  5M  HCI  and  0.6%  sodium  borohydride 
(NaBH,)  in  0.5%  NaOH.  Flow  rates  were  approximately  7  and 
1  mL  min"1  for  samples  and  reagents,  respectively.  Generated 
arsine  was  quantified  at  193.7  run  in  an  air-acetylene  flame 
(Perkin  Elmer  3100).  The  As  detection  limit,  estimated  as 
three  standard  deviations  of  50  absorbance  blank  readings, 
was  3.4  nM.  Free  ion  activities  were  estimated  using  the 
aqueous  chemical  equilibrium  model,  GEOCHEM  (Sposito 
and  Mattigod,  1979)  as  modified  by  Parker  et  al.  (1987). 
Stability  constants  for  pertinent  aqueous  complexes  were  ob- 
tained from  compiled  thermodynamic  equilibrium  data  (Sadiq 
and  Lindsay,  1979;  Sposito  and  Mattigod,  1979;  Nordstrom 
and  May,  1989). 

Each  sample  was  subjected  to  a  sequential  extraction  for 
As  (Johnston  and  Barnard,  1979;  Ganje  and  Rains,  1982) 
based  on  developed  procedures  for  extraction  of  P  (Olsen  and 
Sommers,  1982)  and  organically  bound  metals  (Tessler  et  al., 
1979).  Although  sequential  extractions  yield  only  qualitative 
information  on  As  fractionation,  this  procedure  was  useful  for 
assessing  the  relative  lability  of  As  among  the  different  samples. 
The  method  consisted  of  sequential  extractions  in  the  following 
order:  1  M  NH4C1,  0.1  M  NaOH,  0.3  M  citrate  bicarbonate 
(CB),  0.3  M  citrate  dithionite  bicarbonate  (CDB),  1  M  HCI 
and  30%  H,O2/0.8  M  NH4OAc.  All  samples  were  analyzed 
for  As  with  HG-AAS;  the  CDB  and  H202  extracts  were  also 
analyzed  with  ICP.  The  residual  solids  were  digested  using  a 
four  acid  heat-treatment  (Welsch  et  al..  1990)  and  As  levels 
in  the  resulting  digests  were  determined  using  HG-AAS. 

Selected  grains  of  PT  and  RT  were  affixed  to  Al  stubs  and 
coated  with  C  for  viewing  with  a  scanning  electron  microscope 
(SEM).  The  grains  were  mapped  for  15  elements  using  Energy 
Dispersive  Analysis  of  x-rays  (EDAX)  to  identify  particles 
enriched  in  As,  Fe,  or  S.  Particles  that  contained  high  levels 
of  these  elements  were  analyzed  in  further  detail  to  determine 
semiquantitative  atomic  ratios  using  standardless  analysis  with 
PROZA  correction  (Goldstein  et  al  ,  1992). 


Column  Experiments 

Unsaturated  column  experiments  were  conducted  on  limed 
and  unlimed  PT  and  RT.  Liming  treatments  consisted  of  adding 
a  mix  of  60%  CaCO,/40%  Ca(OH)2  at  concentrations  of  3.0 
and  25.6  g  kg"1  of  PT  and  RT,  respectively.  The  CaCO,/ 
Ca(OH)2  ratio  and  liming  concentrations  were  determined  from 


Table  1 .  Total  arsenic  and  metal  concentrations  of  six  mine  sam- 
ples determined  by  x-ray  fluorescence  (XRF). 


Smelter  stack  deposits 

Reprocessed 

Pond  tailings 
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Element 
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B            C 
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.    _, 

As 

248 

2  358            48 

1  506 

3  421 

2  338 

(  .1 

NDt 

27  900     21  100 

20  700 

14  800 

5  300 

Cu 
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1  331           126 
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2  386 

2  382 

Ee 

13  100 

46  800     42  800 

35  400 

215  800 

1 52  300 

Mn 
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1  283       1  203 
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914 

ND 

Pb 

704 

301       ND 

147 

881 

813 

Zn 

103 

816         205 

466 

341 

671 

t  ND  = 

none  detected. 

acid-base  accounting  as  part  of  a  remedial  investigation  study 
(Schafer  et  al..  1989).  The  FT  (25%  sand,  56%  silt,  19% 
clay)  was  diluted  3:1  with  40  to  200  mesh  SiO;  (Fluka)  to 
assure  adequate  flow,  whereas  the  RT  (95%  sand,  3%  silt, 
2%  clay)  was  left  undiluted.  Columns  (5  cm  dtam.,  20  cm 
long  PVC  pipe)  were  sealed  with  rubber  O-rings  in  polycarbo- 
nate end  caps.  A  plastic  porous  plate  (Soil  Measurement  Sys- 
tems, Tucson,  AZ)  was  inserted  within  each  bottom  end  cap 
to  control  soil  matric  potential.  Soil  columns  were  packed  by 
tapping  lightly  between  successive  soil  additions  and  then  were 
saturated  from  the  bottom  with  5  mM  K2S04.  Upon  complete 
wetting,  a  20  kPa  vacuum  was  applied  to  the  bottom  of  the 
column,  and  5  mM  K2S04  was  delivered  to  the  top  of  the 
column  at  10  mL  h"'  with  a  syringe  pump.  This  flow  rate 
equates  to  a  pore  water  velocity  of  approximately  1.7  cm  h-1. 
The  vacuum  was  adjusted  over  the  course  of  the  experiments 
to  maintain  volumetric  water  contents  (9V)  between  30  to  36% 
and  23  to  28%  for  the  PT  and  RT,  respectively.  These  water 
contents  were  well  below  the  saturated  water  contents  in  each 
respective  column.  Under  these  conditions  the  columns  re- 
mained aerated  and  it  was  assumed  that  soluble  As  was  present 
as  arsenate  (As(V])  rather  than  arsenite  (As[IlI|).  Unlimed 
control  columns  were  sampled  for  approximately  23  and  100 
pore  volumes  for  the  PT  and  RT,  respectively,  then  treated 
with  a  high  pH  influent  (0. 1  M  K.HCO,,  adjusted  to  pH  10  1 
with  KOH,  and  bubbled  with  air  to  maintain  the  pH  level) 
All  effluent  samples  were  analyzed  for  pH,  EC,  As,  and 
alkalinity.  Selected  samples  were  also  analyzed  for  Fe,  Ca, 
S,  Cu,  and  Zn,  and  a  smaller  subset  analyzed  for  K,  Mn,  Al, 
Na,  and  Mg.  Analysis  methods  were  identical  to  those  outlined 
previously.  Free  ion  activities  were  calculated  using 
GEOCHEM  on  those  samples  analyzed  for  the  complete  param 
eter  set,  as  well  as  selected  samples  which  produced  reasonable 
cation/anion  balances  without  data  on  K.  Mn,  Al,  Na,  and 
Mg. 

RESULTS  AND  DISCUSSION 

Sample  Characterization 

Total  As  contents  of  the  six  mine  samples  (Tabic  1 ) 
ranged  from  48  to  3421  mg  kg"1  (10-700-fold  higher 
than  typical  background  levels  of  5  mg  kg"';  Nriagu, 
1994).  Iron  represented  15  to  21%  of  the  RT,  indicative 
of  high  initial  pyrite  (FeS?)  content.  Total  metal  concen- 
trations of  Cu,  Pb,  and  Zn  were  also  generally  elevated, 
consistent  with  the  composition  of  mine  tailings  in  this 
region.  Stabilization  of  these  barren  tailings  through 
revegetation  has  been  a  primary  goal  of  mine  reclamation 
efforts  to  reduce  erosion  and  prevent  migration  of  these 
contaminants  to  receiving  waters. 


JONES  ET  AL 

ARSENIC  TRANSPORT  IN  MINE  TAILINGS 

435 

Table  2.  Soil  water  chemistry 

of  saturated  paste  extracts  for  six 

mine  soils  and  saturation  indicest 

with  respect  to  selected  solid  phases. 

Pond  tailings  (PT) 
A 

Smelter  stack  deposits  (SSD) 

Reprocessed  taili 

ngs  (RT) 

Constituent 

B 

< 

D 

E 

F 

%  Saturation,  mass 

basis 

37 

29 

27 

25 

24 

24 

pH 

4.1 

7.8 

8.0 

7.8 

3.5 

3.9 

E.C.,  dS/m 

1.52 

0.62 

0  35 

0.36 

2.79 

2.56 

DOC,  mMasC 

3  8 

2.5 

17 

0.8 

0.3 

II  4 

Noncarb.  alk.,  mM 

charge 

0 

2.0 

0.3 

0.9 

0 

II 

Cations 

K,  mAf 

0.31 

0.17 

0.11 

0.04 

0.08 

0.19 

Ca,  mAf 

5.19 

2.08 

0.85 

1.14 

13.68 

13.97 

Mg,  mM 

1.28 

1.08 

0.47 

043 

0.41 

040 

Na,  mM 

0.48 

0.80 

0  71 

0.66 

0.10 

0.12 

Al,  uAf 

149 

<7 

<7 

<7 

1067 

360 

Fe,  uAf 

26 

0  4 

<0.30 

<0.30 

375 

K'ri 

Mn,  uAf 

269 

1  8 

<0.05 

<0.05 

71 

)9 

Zn,  uAf 

KM 

0.9 

0.5 

0.2 

74/> 

767 

Cu,  uAf 

510 

:  8 

1  1) 

1.8 

528 

502 

Pb,  uAf 

<0.41 

<0.41 

<0.41 

<0.41 

0.62 

1.20 

Ni,  uAf 

0.97 

<0.37 

<0.37 

<0.37 

6.00 

3.53 

Cr,  uAf 

0  35 

<0.27 

<0.27 

<0.27 

0.56 

0.65 

Cd,  uAf 

1  09 

<0.05 

<0.05 

<0.0S 

0.96 

1.61 

Anions,  mAf 

As 

0.0007 

2.32 

0.02 

0.64 

0.0003 

0.0001 

B 

0.007 

0.005 

0.003 

0.007 

0.018 

0.019 

CI 

0.27 

0.10 

0.14 

0.05 

0.06 

0.05 

HCO, 

0.00 

1.33 

2.11 

0.74 

0.00 

0.00 

NO, 

0.33 

0.40 

0.46 

0.19 

0.07 

0.05 

PO< 

NDJ 

0.03 

ND 

0.01 

ND 

ND 

SO, 

8.26 

1.51 

0.26 

0.67 

18.22 

13.30 

Charge  balance! 

Lcationic  charge 

mAf(  +  ) 

12.36 

5.07 

3.81 

3.40 

20.85 

21.80 

Lanionic  charge, 

mAf(-) 

12.06 

5.13 

2.44 

2.41 

27.19 

21.28 

Net  charge,  mAf 

0.30 

-0.06 

1.37 

0.98 

-6.34 

0.52 

%  error,  as  cationic  chg 

2.4 

-1.1 

35.9 

29.0 

-30.4 

2  4 

Saturation  indices 

AHAsO,) 

-5 

-8 

-10 

-8 

-7 

-  7 

Ca^AsO,); 

-23 

-  5 

-9 

-6 

-26 

-25 

CuKAsO.), 

-9 

-  1 

-  7 

-2 

-13 

-12 

Fe(AsO.) 

-  4 

-  3 

-8 

-5 

-4 

-3 

MgjfAsO.), 

-23 

-4 

-8 

-t, 

-28 

-27 

Mn^AsO.); 

-16 

-4 

-13 

-11 

-22 

-22 

Zn^AsO.fc 

-17 

-7 

-11 

-9 

-20 

-19 

t  Saturation  index  =  log  [IAPtion  activity  product)//^):  K„  values  from  Wagemann  (1977) 

t  ND  =  none  detected. 

§  Charge  balance  determined  from  distribution  of  ion  species  (GEOCHEM). 


The  chemical  content  of  saturated  paste  extracts  for 
the  six  mine  tailing  samples  is  shown  in  Table  2.  The 
SSD  had  saturated  paste  pH  values  ranging  from  7.8  to 

8.0,  whereas  the  PT  and  RT  had  pH  levels  of  3.5  to 

4.1.  As  expected,  soluble  concentrations  of  nearly  all 
metal  ions  were  highest  in  those  samples  with  low  pH, 
with  SO4  being  the  principal  anion.  High  soluble  con- 
centrations of  Zn,  Cu,  Ni,  Cr,  and  Cd  in  the  low  pH 
samples  document  the  need  for  liming  to  reduce  potential 
soluble  metal  migration  and  phytotoxicity.  Conversely, 
the  SSD  had  relatively  low  soluble  metal  concentrations, 
balanced  primarily  with  the  anions  HCO3  ,  SO4  ,  and 
HAsOj ".  Soluble  As  concentrations  for  the  six  soils 
spanned  more  than  four  orders  of  magnitude,  and  were 
positively  correlated  with  pH.  The  relationship  between 
soluble  As  and  pH  among  all  samples  is  consistent  with 
the  pH  dependence  of  As  sorption  reactions  on  oxide 
minerals  and  layer  silicates  (Goldberg,  1986;  Goldberg 
and  Glaubig,  1988).  Total  As  did  not  correlate  well  with 
soluble  As;  for  example,  RT  samples  had  the  highest 
total  As  but  the  lowest  soluble  As. 

Charge  balance  calculations  performed  after 
GEOCHEM  speciation  yielded  errors  ranging  in  magni- 


tude from  1.1  to  35.9%  (Table  2).  A  comparison  of 
ion  activity  products  (IAPs)  for  seven  possible  metal 
arsenates  with  their  respective  literature  Ksp  values 
(Wagemann.  1977),  indicated  that  all  six  saturated  pastes 
were  undersaturated  with  respect  to  metal  arsenates. 
Given  discrepancies  in  literature  Ksp  values  (Robins, 
1981;  Sadiq  and  Lindsay,  1981;  Robins,  1987),  and 
uncertainty  regarding  the  extent  that  true  equilibrium 
was  reached  during  the  24-h  saturated  paste  extraction, 
it  cannot  be  definitively  concluded  from  IAP  calculations 
that  metal  arsenates  did  not  exist  in  these  samples.  Disso- 
lution studies  with  metal  arsenates  (Ca,  Fe,  and  Mg)  show 
that  As  and  metal  activities  can  vary  up  to  approximately 
twofold  after  24  h  from  calculated  equilibrium  activities 
(Robins,  1981).  At  most,  this  could  account  for  errors 
in  log  (IAP/KSp)  of  0.6  to  1.5  for  Ca,  Fe,  Mg  arsenates, 
thus  still  suggesting  the  saturated  pastes  were  undersatu- 
rated with  respect  to  these  metal  arsenates  (Table  2). 
The  saturation  indices  coupled  with  direct  SEM/EDAX 
observations  (discussed  later)  suggest  that  discrete 
metal  arsenate  phases  were  not  controlling  aqueous  As 
activities. 

The  PT  sample  and  the  lower  pH  RT  sample  (Soils 
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Fig.  1.   Percentage  of  total  As  recovered  from  each  sequential  extract- 
ant  for  both  the  reprocessed  tailings  and  pond  tailings. 


A  and  E,  respectively,  Table  2)  were  chosen  for  further 
sample  characterization  and  column  experiments  because 
their  low  pH  has  apparently  prevented  natural  revegeta- 
tion  in  the  field.  Liming  is  therefore  the  likely  remediation 
strategy  (Neuman  et  al.,  1993b).  Sixty  percent  of  As  in 
the  RT  was  found  in  the  residual  fraction  (Fig.  1),  with 
32%  of  the  remainder  in  the  NaOH  fraction  (nonoccluded 
Fe/Al  oxides;  Olsen  and  Sommers,  1982).  Conversely, 
only  9%  of  the  As  in  the  PT  existed  in  the  residual 
fraction,  while  66%  was  found  in  the  NaOH  fraction. 
Although  0. 1  M  NaOH  could  also  extract  As  bound  as 
Fe  or  Al  arsenates,  the  saturated  paste  data  suggested 
that  these  materials  were  undersaturated  with  respect  to 
these  metal  arsenate  phases.  The  large  difference  in 
readily  extractable  As  between  these  two  samples  sug- 
gests a  significant  difference  in  As  lability.  Arsenic  in 
the  PT  sample  may  be  considered  more  susceptible  to 
mobilization  compared  to  the  RT,  as  will  be  discussed 
further. 

Analysis  of  PT  and  RT  using  SEM-EDAX  revealed 
no  solid  phases  concentrated  in  As.  Most  particles  con- 
tained As  near  or  below  the  detection  limit  of  the  instru- 
ment (approximately  0.2%   by  atom),  although  some 
particles  contained  as  much  as  6%  As.  Particles  con- 
taining high  levels  of  Fe,  S,  or  As  were  further  analyzed 
for  possible  identification  of  primary  or  secondary  miner- 
Table  3.  Molar  ratios  of  S/Fe  and  As/Fe  in  particles  of  pond 
tailings  and  reprocessed  tailings  as  determined  by  energy  disper- 
sive analysis  of  x-rays  (EDAX).  For  comparison,  ratios  of  S/ 
Fe  and  As/Fe  are  shown  for  several  possible  mineral  phases. 


Pond 

tailings 

Reprocessi 
S/Fe 

td  tailings 

Particle 

S/Fe 

As/Fe 

As/Fe 

1 

0.05 

0.11 

0.52 

0.01 

2 

0.89 

0.13 

0.25 

0.02 

3 

0.67 

0.04 

0.18 

0.02 

4 

1.80 

0.36 

0.53 

0.02 

5 

0.35 

0.11 

0.16 

0.03 

6 

1.15 

0.08 

0.12 

0.03 

7 

1.04 

0.07 

0.40 

0.04 

8 

0.93 

0.48 

0.31 

0.05 

Mean  (SD) 

0.86(0.53) 

0.17  (0.16) 

0.31  (0.16) 

0.03  (0.01) 

Mineral 

Arsenopyrite 

1 

1 

1 

1 

Pyrite 

2 

- 

2 

- 

Jarosite 

0.67 

- 

0.67 

- 

Scorodite 

- 

1 

- 

1 

als.  The  S/Fe  ratios  in  both  samples  were  generally  well 
below  elemental  ratios  for  pyrite,  and  fell  closer  to 
expected  ratios  for  jarosite  (Table  3).  Although  some  of 
the  PT  particles  had  S/Fe  ratios  near  arsenopyrite,  As/ 
Fe  ratios  were  all  well  below  the  1 : 1  ratio  found  in  either 
arsenopyrite  (FeAsS)  or  scorodite  (FeAs042H20).  It 
is  difficult  to  invoke  from  EDAX  analyses  alone  the 
presence  of  specific  primary  or  secondary  minerals  of 
Fe,  S,  or  As.  The  low  S/Fe  ratios  in  most  of  these 
samples  (as  compared  to  pyrite)  suggest  conservation 
of  Fe  (likely  as  Fe  oxyhydroxides)  relative  to  S.  The 
formation  of  jarosite  would  not  be  unexpected  in  these 
samples,  given  that  it  is  a  common  secondary  mineral 
formed  from  the  oxidation  of  pyrite  (Dixon  and  Weed, 
1977);  however,  we  found  no  evidence  for  the  presence 
of  residual  arsenopyrite  or  secondary  formation  of  scoro- 
dite. It  is  likely,  based  on  the  strong  affinity  of  As  for 
Fe  oxyhydroxides  (Hingston  et  al.,  1971;  Goldberg, 
1986;  Fuller  et  al.,  1993),  that  the  majority  of  As  was 
sorbed  to  these  minerals. 

Column  Experiments 

Effluent  As  concentrations  of  the  unlimed  PT  column 
reached  steady-state  levels  near  10"65  M  within  several 
pore  volumes  (Fig.  2A).  Liming  the  PT  resulted  in  an 
effluent  As  concentration  that  was  approximately  two 
orders  of  magnitude  higher  than  the  unlimed  column. 
After  13  pore  volumes,  effluent  from  the  limed  PT  col- 
umn peaked  at  approximately  10"4  M\  400-fold  higher 
than  in  the  unlimed  column.  The  amount  of  As  mobilized 
from  the  limed  column  during  the  experiment  (26  pore 
volumes)  represented  28%  of  the  initial  column  As. 
Effluent  pH  values  of  the  unlimed  PT  column  were 
relatively  stable  at  approximately  4.5  to  5  after  several 
pore  volumes.  Liming  resulted  in  an  initial  effluent  pH 
near  8,  followed  by  a  gradual  climb  to  pH  values  near 
10.0.  Increases  in  effluent  pH  of  the  limed  column  corre- 
sponded to  increases  in  effluent  As  (Fig.  2A),  consistent 
with  the  pH  dependence  of  As  sorption  reactions  on 
oxide  minerals  (Goldberg,  1986).  Addition  of  0.1  M 
KHCO3  (pH  10. 1)  to  the  unlimed  PT  column  at  23  pore 
volumes  resulted  in  an  increase  in  effluent  pH  and  As; 
however,  after  three  additional  pore  volumes,  the  hydrau- 
lic conductivity  decreased  to  near  0,  likely  due  to  disper- 
sion of  colloidal  material,  and  no  further  samples  could 
be  obtained. 

Effluent  As  concentrations  remained  near  10  78  M  at 
pH  levels  near  3.5  in  the  unlimed  RT  column  for  100 
pore  volumes  (Fig.  2B).  The  liming  treatment  resulted 
in  an  increase  in  effluent  pH  to  values  near  8  and  effluent 
As  to  values  near  10~7  M.  Although  the  liming  amend- 
ment increased  effluent  As  concentrations  by  almost  one 
order  of  magnitude  in  the  RT  sample,  As  concentrations 
ranged  from  two  to  three  orders  of  magnitude  lower 
than  those  found  in  the  PT  after  liming.  The  total  As 
content  of  the  PT  column  was  only  62  mg  kg  '  (following 
dilution  with  SiOj),  compared  with  approximately  3400 
mg  kg "'  for  the  RT  column.  Therefore,  the  soluble/total 
As  ratio  was  approximately  5000-fold  higher  in  the  limed 
PT  column  than  in  the  RT  column. 
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Fig.  2.  Arsenic  concentrations  and  pH  levels  in  column  effluent  for  limed  and  unlimed  treatments  of  pond  tailings  (Fig.  2A)  and  reprocessed 
tailings  (Fig.  2B).  The  arrow  indicates  the  point  of  addition  of  0.1  M  KHCOj  (pH  10.1)  influent  to  the  unlimed  column. 


Differences  in  effluent  As  concentrations  between  the 
two  limed  samples  were  consistent  with  differences  in 
the  "lability"  of  bound  As  (Fig.  1).  The  NaOH-extractable 
As  fraction  is  thought  to  represent  As  associated  with 
Al  or  Fe  and  would  likely  be  more  susceptible  to  desorp- 
tion  with  increases  in  pH  following  liming.  Presumably, 
the  large  residual  fraction  in  the  RT  (60%)  would  not 
be  as  "labile"  after  liming.  Differences  in  As  mobilization 
between  the  two  samples  may  have  also  been  influenced 
by  differences  in  pH  values  after  liming.  Effluent  pH 
values  in  the  PT  column  reached  9.8  within  7  pore 
volumes  and  remained  near  9  for  almost  20  pore  volumes. 
Conversely,  pH  values  in  the  RT  column  remained  near 
8  for  the  majority  of  the  column  experiment.  Based  on 
the  pH  dependence  of  As  sorption  reactions  on  oxide 
minerals  (Goldberg,  1986),  higher  pH  values  in  the  PT 
column  likely  contributed  to  more  As  desorption.  In  fact, 
when  the  pH  of  the  unlimed  RT  column  was  raised  to 
near  10,  effluent  As  increased  from  approximately  10~7 
to  10"3  2  M  within  20  pore  volumes.  The  fact  that  a  high 


ionic  strength/high  pH  influent  solution  (0. 1  M  KHCOj) 
was  necessary  to  mobilize  substantial  amounts  of  As 
from  the  RT  column  is  consistent  with  the  higher  residual 
As  fraction  in  this  sample.  Results  from  the  column 
experiments  presented  here  are  consistent  with  the  pH 
dependence  of  surface  complexation  reactions  of  As  on 
solid  phase  minerals  such  as  goethite,  ferrihydrite,  gibb- 
site,  Al  hydroxides,  and  layer  silicates. 

Prior  to  liming,  it  was  demonstrated  that  IAPs  of  PT 
and  RT  saturated  paste  extracts  were  undersaturated  with 
respect  to  several  metal  arsenate  solid  phases  (Table  2); 
however,  after  liming,  changes  in  pore  water  chemistry 
necessitated  a  reanalysis  of  IAPs,  especially  considering 
the  potential  for  Ca-arsenate  to  form  in  the  presence  of 
calcite  (Sadiq  et  al.,  1983).  Solubility  curves  for  both 
scorodite  and  Ca-arsenate  were  compared  to  selected 
effluent  data  for  each  of  the  three  treatments  (Fig.  3). 
All  effluent  solutions  were  found  to  be  undersaturated 
with  respect  to  these  two  solid  phases.  Although  the 
possibility  of  an  Fe-arsenate  phase  with  lower  solubility 
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Fig.  3.  Solubility  diagrams  for  scorodite  (FeAs04  2H2O)  and 
Caj(AsO<)2  compared  to  activities  of  AsOJr,  Fe3*,  and  Ca!+  in 
column  effluent.  Data  points  below  the  solubility  lines  indicate  these 
solutions  are  undersaturated  with  respect  to  the  given  solid  phase. 


than  that  assumed  for  scorodite  cannot  be  ruled  out 
(Nordstrom  and  Parks,  1987),  such  phases  have  not  been 
commonly  observed  in  other  studies.  In  addition,  other 
possible  solid  phase  Ca-arsenates  such  as  Ca(HiAs04)2, 
CaHAsCX,  and  Ca5H2(As04)4  can  exist  in  a  Ca-As-H:0 
system.  Each  of  these  solid  phases  is  more  soluble  than 
Ca3(As04)2  (Nishimura  et  al.,  1987),  hence  the  effluent 
solutions  were  undersaturated  with  respect  to  these  solid 
phases  as  well.  These  data,  combined  with  the  EDAX 
results,  support  the  contention  that  the  behavior  of  As 
after  liming  is  not  controlled  by  equilibrium  with  respect 
to  solid  phase  arsenates,  but  is  related  to  the  pH  depen- 
dence of  As  sorption  reactions  typically  observed  on 
oxide  minerals. 


CONCLUSIONS 

Characterization  of  several  representative  mine  tailing 
samples  showed  that  soluble  As  levels  were  not  correlated 
with  total  As,  but  were  more  strongly  correlated  with 
soil  solution  pH.  Furthermore,  saturated  pastes  of  all 
mine  tailing  samples  were  undersaturated  with  respect 
to  known  metal  arsenate  solid  phases.  The  mobility  of 
As  in  two  selected  mine  tailing  samples  was  shown  to 
increase  substantially  in  unsaturated  laboratory  columns 
following  a  recommended  liming  treatment  (CaCOj/ 
Ca[OH]:).  Adsorption-desorption  reactions  on  oxide 
minerals  are  believed  to  be  responsible  for  the  pH  depen- 
dence of  soluble  As  concentrations,  rather  than  precipita- 
tion-dissolution reactions  of  solid  metal  arsenates.  The 
concentration  of  As  in  column  effluent  after  liming  was 
not  correlated  with  total  As  in  the  two  samples  chosen 
for  detailed  column  experiments.  Sequential  extractions 
showed  substantial  differences  in  "labile"  As  between 
the  two  samples,  consistent  with  results  from  column 
experiments  after  liming.  Consequently,  the  distribution 


of  soil-bound  As  is  important  for  determining  the  suscep- 
tibility of  mine  tailings  to  As  mobilization  following 
liming.  In  addition,  the  pH  of  the  post-limed  tailings  is 
important  for  determining  the  amount  of  As  mobilized 
after  liming.  Remediation  of  acidic  mine  tailings  or  other 
As-contaminated  soils  using  lime  amendments  should  be 
evaluated  with  respect  to  the  potential  effects  on  As 
mobilization,  especially  at  contaminated  sites  hydrau- 
lically  connected  to  surface  or  groundwaters. 
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ABSTRACT 

The  effects  of  pH  and  phosphate  (PO,)  competition  on  arsenate 
(AsO,)  sorption  by  soils  and  aquifer  materials  have  been  previously 
investigated  in  batch  studies,  however,  little  has  been  reported  on 
processes  that  influence  AsO,  transport  through  porous  media.  Conse- 
quently, the  objective  of  this  study  was  to  determine  the  effects  of 
pH  and  PO,  on  the  transport  of  AsO,  under  conditions  where  reaction 
and  mass  transfer  rates  may  be  controlling  pore-water  AsO,  concen- 
trations. Saturated  column  transport  experiments  (pore  water 
velocity  =  1  cm  h  ')  were  performed  using  a  sand  containing/rev  Fe 
oxides  in  which  a  pulse  of  "AsO,  ( 133  (i  M )  was  applied  in  the  presence 
of  0,  13.4,  134,  and  1340  p.A/  PO,  at  pH  values  of  4.5,  6.5,  and  8.5. 
Breakthrough  curves  (BTCs)  of  AsO,  exhibited  low  recovery  and 
significant  tailing  at  pH  values  of  4.5  and  6.5  in  the  absence  of  PO,. 
Competition  between  P04  and  AsO,  for  sorption  sites  resulted  in 
increased  AsO,  mobility;  however,  even  in  the  presence  of  PO,  levels 
representing  more  than  100%  of  maximum  sorption  capacity  for  the 
columns,  significant  amounts  of  AsO,  remained  sorbed  to  the  sand. 
Application  of  a  continuous  PO,  pulse  to  a  column  containing  sorbed 
AsO,  resulted  in  an  increase  in  AsO,  recovery  (>35%);  however, 
even  after  total  PO,  loading  exceeded  the  column  capacity  by  more 
than  two  orders  of  magnitude,  approximately  40%  of  the  applied 
AsO,  remained  sorbed  to  the  sand  indicating  that  rates  of  AsO, 
desorption  play  an  important  role  in  transport  of  AsO,  through  po- 
rous media. 


The  pH  dependence  of  As04  adsorption  has  been 
investigated  on  numerous  solid  phases  including 
amorphous  Al  hydroxide  (Anderson  et  al.,  1976),  amor- 
phous Fe  hydroxide  (Pierce  and  Moore,  1982),  goethite 
and  gibbsite  (Goldberg,  1986;  Hingston  et  al.,  1971), 
kaolinite,  montmorillonite,  calcite,  and  quartz  (Gold- 
berg and  Glaubig,  1988;  Frost  and  Griffin,  1977;  Xu  et 
al.,  1991).  It  is  well  established  that  adsorption  of  As04 
on  Fe  and  Al  hydroxides  (including  goethite  and  gibbs- 
ite) increases  with  decreasing  pH  to  pH  values  near  3. 
Arsenate  sorption  on  calcite  and  layer  silicate  minerals 
generally  shows  a  maximum  in  the  pH  range  of  4  to  6. 
Adsorption  of  As04  onto  clean,  crystalline  quartz  was 
demonstrated  to  be  negligible  at  pH  values  above  3  (Xu 
et  al.,  1991).  The  existing  adsorption  data  show  that  in 
natural  systems,  As04  is  preferentially  adsorbed  to  Fe 
and  Al  oxide  minerals  compared  to  kaolinite,  montmo- 
rillonite, calcite,  or  quartz  (Fordham  and  Norrish,  1979; 
Jacobs  et  al.  1970;  Johnston  and  Barnard,  1979;  Livesey 
and  Huang,  1981). 

Research  has  also  shown  substantial  competition  be- 
tween AsO,  and  P04  sorption  on  a  variety  of  soils  and 
oxide  minerals  (Goldberg,  1986;  Hingston  et  al.,  1971; 
Peryea,  1991;  Roy  et  al.,  1986a,b).  In  soils  that  have 
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been  previously  contaminated  with  lead  As04  pesti- 
cides, application  of  P04  fertilizers  has  been  shown  to 
result  in  increased  As  mobility  (Peryea,  1991).  Hingston 
et  al.  (1971)  modeled  the  competition  between  AsO, 
and  PO„  on  both  goethite  and  gibbsite  by  a  modified 
Langmuir  isotherm  using  selectivity  coefficients  derived 
from  experimental  data  involving  binary  mixtures  over 
a  wide  range  of  concentrations.  This  approach  has  been 
extended  to  natural  soils  as  well  (Roy  et  al.,  1986a,b). 
Goldberg  (1986)  successfully  used  a  constant  capaci- 
tance model  to  predict  the  competition  between  AsO, 
and  P04  on  gibbsite  over  a  range  of  As04  and  P04 
solution  concentrations,  with  surface  complexation  con- 
stants determined  from  a  binary  system  in  which  AsO, 
and  P04  were  present  in  equimolar  amounts. 

While  there  has  been  substantial  research  on  the  ef- 
fects of  pH  and  P04  competition  on  As04  adsorption 
in  batch  studies,  there  is  little  data  regarding  the  effects 
of  these  processes  on  transport  of  As04  through  porous 
media.  Furthermore,  while  these  processes  have  been 
quantified  and  modeled  in  batch  systems,  such  systems 
typically  are  well  mixed  with  high  solution  to  solid  ratios, 
and  represent  near  equilibrium  conditions.  It  is  unclear 
the  extent  to  which  such  data  may  be  extrapolated  to 
transport  environments  where  reaction  and  mass  trans- 
fer rates  may  be  controlling  pore  water  As04  concentra- 
tions (Darland  and  Inskeep,  1997).  Consequently,  the 
objectives  of  this  research  were  to  determine  the  (i) 
effect  of  pH  on  As04  transport  through  a  representative 
sand,  (ii)  extent  to  which  simultaneous  competition  be- 
tween As04  and  PG4  for  adsorption  sites  affects  As04 
transport,  (iii)  effect  of  pH  on  the  competition  between 
P04  and  As04  for  adsorption  sites  during  transport,  and 
(iv)  extent  to  which  modification  of  pore  water  P04 
concentration  could  increase  As04  mobility  in  contami- 
nated aquifer  sands.  The  focus  of  this  study  was  to  eluci- 
date processes  important  in  the  transport  of  As  under 
oxidizing  conditions  where  As(V)  is  the  predominant 
species.  Appropriate  mathematical  expressions  neces- 
sary for  modeling  the  transport  of  As  under  column 
conditions  are  discussed  in  a  separate  study  (Darland 
and  Inskeep,  1997). 

MATERIALS  AND  METHODS 

Reagents  and  Analysis 

Sand  of  mixed  mineralogy  (quartz,  orthoclase,  and  plagio- 
clase  feldspars)  was  obtained  from  the  Unimin  Corporation 
(Emmett,  ID),  sieved  to  a  particle-size  range  of  0.25  to  0.50 
mm,  acid  washed  three  times  with  0.05  M  HC1,  then  rinsed 
with  double  deionized  water  (DD-H40)  until  the  conductivity 
of  the  rinse  had  stabilized  at  <2  u.S.  The  sand  was  then  oven 
dried  at  105°C,  and  characterized  for  surface  area  and  chemical 


Abbreviations:  BTCs,  breakthrough  curves;  LS,  liquid  scintillation; 
IC,  ion  chromatography;  AAM,  ascorbic  acid  method;  SA,  surface 
area;  CD,  cilrate-dithionite;  EXAFS,  Extended  X-ray  Absorption 
Fine  Structure  Spectroscopy;  lAPs,  ion  activity  products. 
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composition  (Tabic  1).  In  addition,  an  acid  purified  quartz 
sand  (SiO,)  was  obtained  (Fluka  Chemical  Co.,  Ronkonkoma, 
NY)  presieved  to  a  particle-size  distribution  of  0.074  to  0.42 
mm  and  used  in  selected  As04  transport  experiments.  All 
chemical  impurities  in  the  purified  quartz  sand  were  below 
detection,  as  given  by  the  manufacturer. 

All  As04  solutions  were  prepared  with  the  same  source  of 
analytical  grade  Na;HAs(J4  salt,  and  diluted  with  DD-HjO. 
Arsenate  concentrations  were  measured  by  three  different 
methods  depending  on  the  experiments  performed.  For  the 
batch  studies,  As04  was  measured  by  atomic  absorption  with 
continuous  hydride  generation  (HG-AAS).  For  the  transport 
studies,  radiolabeled  ''As04  (Los  Alamos  National  Labora- 
tory, Los  Alamos,  NM)  was  used  in  conjunction  with  liquid 
scintillation  (LS)  counting.  Ion  chromatography  (IC)  was  used 
(Column  =  Dionex  AS4A-SC;  Eluant  =  2.5  mM  Na2CO,/2.66 
mM  NaHCO,;  Flowrate  =  2  mL/min)  periodically  to  verify 
concentrations  determined  by  LS,  and  to  verify  that  no  reduc- 
tion of  arsenate  to  arsenite  had  occurred. 

All  P04  solutions  were  prepared  with  the  same  source  of 
analytical  grade  K2HP04  salt,  and  diluted  with  DD-HjO.  Phos- 
phate concentrations  were  measured  by  two  different  methods 
depending  on  the  experiments  performed.  For  batch  sorption 
studies.  P04  was  measured  by  the  standard  ascorbic  acid 
method  (AAM)  (APHA,  1989).  For  transport  studies,  As04 
interference  using  the  AAM  necessitated  the  use  of  IC  (Dio- 
nex AS4A-SC  column)  for  P04  analysis. 


adsorbed  concentration  to  the  concentration  in  solution)  was 
not  reached  in  both  cases  until  at  least  96  h. 

Batch  Isotherms 

Batch  adsorption  isotherms  for  As04  onto  the  sand  were 
performed  at  three  initial  P04  concentrations  (0, 134,  and  1340 
P-M)  in  35-mL  polypropylene  centrifuge  tubes,  using  a  1:1 
solid/solution  ratio  (IS  g  solid:15  mL  solution).  Initial  concen- 
trations of  0,  1.3,  6.7,  13,  67.  133,  and  266  p.M  As04  were  used 
at  each  initial  P04  concentration.  Samples  were  shaken  for 
96  h,  and  subsequently  decanted  and  filtered  through  0.45  u.m 
nylon  filters  prior  to  analysis  for  As04  by  HG-AAS.  Analysis 
of  selected  samples  using  IC  indicated  that,  within  the  limits 
of  detection,  all  of  the  total  soluble  As  was  present  as  As04 
species  (i.e.,  As(V)).  Adsorbed  concentrations  for  As04  were 
determined  by  difference  between  initial  and  final  measured 
solute  concentrations.  Adsorption  isotherms  at  each  initial 
P04  concentration  were  fit  to  the  Langmuir  model  using  non- 
linear regression  according  to  the  equation: 


r„ 


KL    C 

1  +  KL  ■  C 


where  S  is  adsorbed  concentration  (u.mol  kg1),  C  is  solution 
concentration  (u.mol  L-'),  KL  is  an  empirical  affinity  parame- 
ter (L  u.mol-1),  and  T,,  is  maximum  adsorption  density 
(u.mol  kg-1). 


Kinetic  Adsorption  Study 

Separate  batch  studies  were  performed  to  evaluate  the  ad- 
sorption kinetics  of  As04  and  P04  onto  the  sand,  and  to  deter- 
mine an  equilibrium  time  required  for  batch  adsorption  iso- 
therm experiments.  The  batch  kinetic  adsorption  experiments 
were  performed  in  15-mL  polyethylene  screw  top  bottles,  and 
consisted  of  a  1:1  solid/solution  ratio  (5  gsolid:5  mL  solution), 
using  initial  concentrations  of  150  u.A/  As04  and  360  u.Af  P04. 
All  solutions  were  prepared  in  a  background  of  0.01  M  KG. 
Samples  were  collected  at  time  intervals  of  2,  4,  8,  24,  48,  96, 
and  168  h,  and  were  decanted  and  filtered  through  0.45  u.m 
nylon  filters  prior  to  analysis  (As04  by  HG-AAS,  P04  by 
AAM).  Adsorbed  concentrations  were  determined  by  differ- 
ence between  initial  and  final  measured  solute  concentrations. 
Approximately  97%  of  the  equilibrium  adsorbed  As04  con- 
centration was  reached  within  24  h.  The  P04  adsorption  had 
only  reached  92%  of  the  equilibrium  concentration  within 
24  h    Equilibrium  (as  defined  by  constancy  of  the  ratio  of 

Table  1.   Physical  and  chemical  characteristics  of  sand  used  in 
column  transport  experiments. 


Element 


Total  content 


gli*  '1 

Si 

368 

K 

in 

\l 

38 

(  a 

5.7 

N.i 

3.1 

l> 

0.54 

Mn 

•  0.02 

P 

0.02 

Fe 

4*Free"§ 

0.10 

Al-' 

"Free" 

(I.I  II 

Par 

tide  size: 

0.25-050  mm 

Surf 

face  area:*il 

390 

♦  0.03  nv  kg   ' 

CEOS 

2.7 

x   10 

1  mot,  charge  kg 

stdev 

(5)1 

(4) 

(3) 

(0.4) 
(0.9) 
(0.25) 


((1.(11) 


t  Total  chemical  composition  t>y  Li-mctaboratc  fusion. 

t  Deviations  represent  standard  deviations  of  three  replicates. 

§  Free  Fe  and  Al  oxides  determined  by  cilrale-dilhionite  extraction. 

*j]  Surface  area  determined  using  triple-point  BET-Nj  isotherms. 

#  CEC  =  cation  exchange  capacity. 


Column  Transport  Experiments 

Column  transport  experiments  were  performed  on  the  sand 
under  saturated  conditions  at  three  different  pH  values  (4.5, 
6.5,  and  8.5)  for  a  single  As04  concentration  (133  u-M)  in 
conjunction  with  four  different  P04  concentrations  represent- 
ing molar  P/As  ratios  of  0,  0.1,  1.0.  and  10.  Transport  experi- 
ments were  performed  using  polycarbonate  columns  (2.8  cm 
diam,  8  cm  length)  containing  76.5  g  sand,  resulting  in  a  bulk 
density  of  1.55  g  cm"'  and  a  calculated  porosity  of  0.41.  The 
end  caps  of  the  column  consisted  of  fine  nylon  mesh  (<0.1 
mm)  supported  by  a  Teflon  plate  drilled  with  1-mm  holes  on 
3-mm  centers.  A  syringe  pump  (Soil  Measurement  Systems. 
Tucson,  AZ)  was  used  to  provide  flow  to  the  columns. 

Typical  transport  experiments  were  performed  as  follows. 
After  initial  saturation  and  conditioning  with  0.01  M  KG  to 
the  desired  pH,  approximately  one  pore  volume  of  H:0  (~500 
Bq  mL-1)  in  0.01  M  KG  was  pulsed  through  the  columns 
(direction  of  flow  =  bottom  to  top)  at  an  average  pore  water 
velocity  of  1  cm  h-1  (8  h  per  pore  volume;  0.25  mL  pump 
volumes  at  a  frequency  of  350  s),  followed  by  0.01  M  KG. 
After  complete  recovery  of  the  'H:0  pulse  (~6  pore  volumes 
total  flow),  one  pore  volume  of  a  combined  As04/P04  solution 
(in  0.01  M  KG)  was  applied,  followed  by  0.01  M  KG.  Four 
separate  runs  were  performed  at  each  pH  value,  keeping  the 
As04  concentration  constant  ( 133  u.M;  M650  Bq  mL  '  ;'As), 
and  varying  the  concentration  of  POj  (0.  13.4,  134,  and  1340 
U.A/)  The  experiments  were  continued  for  approximately  20 
pore  volumes  (approximately  7  d)  until  the  tailing  of  As04  in 
columns  effluent  was  well  established  and  As04  concenlra 
tions  were  near  baseline  levels 

The  native  pH  of  column  effluent  using  unbuffered  0.01 
M  KG  was  approximately  4.5.  A  carbonate-buffered  eluant 
solution  (pCO,  =  0.0003  atm)  in  0.01  M  KG  was  used  to 
obtain  constant  column  effluent  pH  values  of  6.5  and  8.5. 
The  buffered  0.01  M  KG  solutions  were  prepared  by  pH 
adjustment  with  KOH  under  conditions  of  continuous  sparg- 
ing with  air  until  the  pH  of  the  eluant  solutions  had  stabilized. 
A  separate  column  experiment  was  performed  at  pH  =  8.5 
under  N:  (g)  to  determine  if  the  HCO,/CO\-  concentrations 
in  the  CO,(g)  buffered  pH  8.5  systems  significantly  affected 
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the  transport  of  AsO„  via  HC03"  or  CO5  competition.  Efflu- 
ent pH  values  were  measured  with  an  in-line  pH  electrode  to 
prevent  sample  exposure  to  atmospheric  C02  prior  to  pH 
measurement.  The  experimental  protocol  was  identical  to  the 
previous  columns,  and  used  the  same  As04  concentration  (133 
u,M;  ~1650  Bq  mL"1  "As)  in  the  absence  of  P04. 

An  additional  column  experiment  was  performed  to  assess 
the  effect  of  P04  on  As04  desorption  in  a  column  already 
loaded  with  an  As04  pulse.  The  column  was  conditioned  (0.01 
M  KCI)  and  loaded  with  a  one  pore  volume  pulse  of  143  (jl/W 
As04  (~1650  Bq  mL"1  "As)  at  a  pH  of  4.5.  After  20  pore 
volumes,  the  0.01  M  KCI  eluant  was  replaced  with  an  eluant 
containing  1420  u-M  P04  (10:1  molar  ratio  of  P/As)  in  0.01  M 
KCI.  The  experiment  was  continued  for  an  additional  15  pore 
volumes  until  sufficient  As  tailing  had  been  established. 

For  comparison  to  the  reactive  sand  used  above,  a  transport 
study  was  also  performed  using  acid  purified  quartz  sand  at 
an  unbuffered  pH  of  approximately  7.0.  The  column  was  con- 
ditioned with  0.01  M  KCI,  and  loaded  with  a  1  pore  volume 
pulse  of  a  133  \lM  As04  solution  and  eluted  according  to  the 
methods  discussed  previously.  The  experiment  was  continued 
until  sufficient  tailing  had  been  established  (~10  pore  vol- 
umes), although  As04  recovery  of  >95%  was  observed  in 
column  effluent  within  five  pore  volumes. 

All  column  effluent  was  passed  through  an  in-line  0.45  jj.ni 
nylon  filter,  then  collected  in  glass  test  tubes  using  a  fraction 
collector  (approx.  eight  samples  per  pore  volume).  Effluent 
samples  were  analyzed  for  pH,  followed  by  "AsO,  using  LS. 
Ion  chromatography  was  used  to  analyze  P04  concentrations 
and  to  verify  reduction  of  arsenate  toarsenite  had  not  occurred 
during  transport.  After  termination,  the  sand  from  the  col- 
umns was  extracted  with  a  mixture  of  6  M  HC1  and  3  M  HNO, 
to  desorb  the  remaining  adsorbed  As04,  and  the  resulting 
extract  analyzed  for  "As  using  LS  to  determine  mass  balance 
of  applied  "AsOj.  In  all  cases,  the  total  recovery  of  applied 
"As04  was  >95%. 

Determination  of  Column  Transport  Parameters 

Individual  sample  volumes  from  the  transport  experiments 
were  determined  by  the  mass  of  the  collected  samples,  and 
the  water  flux  was  determined  as  the  average  value  of  sample 
volume  divided  by  the  collection  interval.  Pore  water  velocity 
was  then  determined  as  the  water  flux  divided  by  volumetnc 
water  content  (q  6V  "'),  as  determined  by  the  equation  8V  = 
1-  (p/2.65)  assuming  a  particle  density  of  2.65  g  cm ",  and 
complete  saturation  in  the  columns,  where  p  is  the  bulk  density 
(g  cm"3).  Relative  concentration  (C/C0)  was  determined  as 
the  ratio  of  the  collected  3H,0  or  As04  concentration  to  the 
applied  'H20  or  As04  concentration  in  the  pulse,  and  the 
number  of  pore  volumes  was  determined  as  the  ratio  of  cumu- 
lative collected  volume  to  the  liquid  volume  of  a  saturated 
column.  Breakthrough  curves  for  'H:0  were  fit  to  the  one- 
dimensional  convection  dispersion  equation,  with  the  com- 
puter program  CXTFIT  (Parker  and  van  Genuchten,  1984), 
and  were  used  to  verify  that  no  transport  related  nonequilib- 
rium  existed  in  any  of  the  column  runs.  Dimensionless  BTCs 
for  As04  were  plotted  using  relative  concentration  (CIC„)  vs. 
pore  volumes  of  the  column  outflow.  Cumulative  recovery  of 
applied  AsO„  collected  in  the  effluent  is  also  presented  on 
all  BTCs. 

RESULTS  AND  DISCUSSION 

Batch  Adsorption 

The  kinetics  of  AsO„  and  P04  adsorption  onto  the 
sand  showed  that  at  least  96  h  was  needed  to  reach 
equilibrium  for  both  As04  and  P04  adsorption  in  a 


well-mixed  batch  system  (Fig.  1).  This  is  consistent  with 
results  from  other  studies  showing  a  slow  approach  to 
equilibrium  for  As04  adsorption  onto  amorphous  Al 
hydroxides  (Anderson  et  al.,  1976)  and  ferrihydrite  (Ful- 
ler et  al.,  1993).  Similarly,  P04  adsorption  onto  Fe  and  Al 
hydroxides  has  been  shown  to  exhibit  a  slow  approach  to 
equilibrium  (Bolan  et  al.,  1985).  Slow  adsorption  kinet- 
ics have  been  previously  attributed  to  diffusion  limita- 
tions within  pore  spaces  of  small  aggregated  particles 
in  the  case  of  P04  adsorption  onto  ferrihydrite  (Willett 
et  al.,  1988).  The  actual  chemical  step  in  the  adsorption 
process  of  As04  and  P04  onto  Fe  and  Al  oxides  and 
clay  minerals  is  quite  rapid,  generally  on  the  order  of 
milliseconds  (Grossl  and  Sparks,  1995;  Amacher,  1991; 
Sparks  and  Zhang,  1991).  The  slow  approach  to  equilib- 
rium observed  in  this  experiment  suggests  that  the  kinet- 
ics of  As04  and  P04  adsorption  onto  the  sand  were 
limited  by  mass  transfer,  such  as  film  and/or  interparticle 
diffusion.  The  sand  used  in  the  current  study  had  a  N2(g) 
surface  area  (SA)  of  0.39  m2  g"1  (Table  1),  which  is 
more  than  50  times  the  estimated  geometric  surface 
area  (spherical).  The  higher  N2(g)-SA  indicates  a  large 
degree  of  surface  roughness  ( Anbeek,  1992),  suggesting 
that  micropores  existed  on  the  sand  particles. 

Adsorption  isotherms  for  As04  onto  the  sand  in  the 
absence  and  presence  of  P04  exhibited  distinct  Langmu- 
irian  behavior  (Fig.  2).  The  observed  data  were  modeled 
by  best  fit  Langmuir  isotherms  (nonlinear  regression), 
yielding  maximum  site  densities  (IYs)  for  As04  of  168, 
112,  and  7  u.mol  kg _1  in  the  presence  of  0, 134,  and  1340 
p.M  initial  P04  (0,  1:1,  10:1  molar  ratios),  respectively. 
The  presence  of  a  1:1  molar  ratio  of  P/As  resulted  in  a 
decrease  in  the  maximum  As04  adsorption  density  of 
33%  compared  to  the  0  \x.M  P  treatment.  In  the  presence 
of  a  10:1  molar  P/As  ratio,  maximum  As04  adsorption 
was  reduced  by  96%.  The  influence  of  P04  competition 
on  As04  adsorption  by  the  sand  (Fig.  2)  is  consistent 
with  As04/P04  competition  reported  in  other  studies 
(Hingston  et  al.,  1971;  Goldberg,  1986). 

Previous  studies  on  the  adsorption  of  As04  by  Fe 
hydroxides  have  shown  that  the  maximum  adsorption 
density  of  As04  onto  ferrihydrite  is  an  order  of  magni- 
tude greater  than  onto  goethite  (Pierce  and  Moore, 
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Fig.  2.  Adsorption  isotherms  for  arsenate  (As04)  in  the  absence  and 
presence  of  134  nA/  and  1340  (iA/  phosphate  (P04)  at  a  1:1  solid/ 
solution  ratio  and  a  pll  range  of  4  to  4.6  (all  in  0.01  M  KCI).  Solid 
lines  represent  best-fit  Langmuir  isotherms,  where  T0  =  maximum 
adsorption  density  and  A',   =  empirical  affinity  parameter. 

1982;  Grossl  and  Sparks,  1995,  Hingston  et  al.,  1971). 
At  a  pH  of  6,  Pierce  and  Moore  (1982)  reported  a 
maximum  site  density  for  ferrihydrite  of  0.12  mol  As04/ 
mol  Fe,  although  no  surface  area  was  reported  for  the 
ferrihydrite  samples  used.  At  the  same  pH,  Grossl  and 
Sparks  (1995)  reported  a  maximum  site  density  for  goe- 
thite  of  0.0084  mol  As04/mol  Fe  with  a  sample  surface 
area  of  50  m2/g,  while  Hingston  et  al.  (1971)  reported 
a  maximum  site  density  for  goethite  of  0.012  mol  AsOV 
mol  Fe  with  a  sample  surface  area  of  60  nr/g. 

Results  of  the  citrate-dithionite  (CD)  extraction  indi- 
cated that  the  sand  used  in  the  current  study  contained 
0.01%  by  weight  free  Fe,  presumably  as  Fe  oxides  and 
hydroxides  (Table  1 ).  If  it  is  assumed  that  As04  adsorp- 
tion by  the  sand  was  a  result  of  coordination  with  surface 
sites  on  Fe  oxides,  a  reasonable  assumption  given  the 
strong  preference  of  Fe  oxides  for  As04  compared  to 
other  soil  phases  (Fordham  and  Norrish,  1979;  Jacobs 
et  al.  1970;  Johnston  and  Barnard,  1979;  Livesey  and 
Huang,  1981),  then  the  maximum  site  density  at  pH  4.5 
for  the  sand  used  in  the  current  study  was  0.095  mol 
AsOVmol  free  Fe.  Pierce  and  Moore  (1982)  reported 
site  densities  on  ferrihydrite  of  0.16  and  0.12  mol  AsOV 
mol  Fe  at  pH  values  of  4  and  5,  respectively,  which 
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is  slightly  greater  than  the  calculated  values  for  this 
experiment.  However,  Pierce  and  Moore  (1982),  used 
a  batch-produced  ferrihydrite  in  a  well-mixed  system, 
which  would  likely  have  a  higher  surface  area  than  the 
free  Fe  oxides  present  on  the  sand  grains  used  in  the 
current  study. 

Effects  of  pH  on  Arsenic  Transport 

The  transport  of  As04  at  pH  values  of  4.5  and  6.5 
exhibited  significant  retardation,  tailing  and  poor  recov- 
ery of  applied  As04  (Fig.  3).  In  contrast,  the  transport 
of  AsG"4  at  pH  8.5  was  more  rapid  and  resulted  in  a 
more  symmetrical  BTC  (Fig.  3).  After  10  pore  volumes, 
the  percent  of  applied  As04  recovered  in  the  column 
effluent  was  35.6,  39.2,  and  100  for  the  columns  at  pH 
4.5,  6.5,  and  8.5,  respectively.  None  of  the  predominant 
mineral  phases  present  in  the  sand  (quartz,  feldspars) 
would  be  expected  to  sorb  appreciable  quantities  of 
As04.  For  comparison,  the  transport  of  As04  through 
acid  purified  silica  sand  exhibited  essentially  no  retarda- 
tion or  tailing  with  As04  recovery  in  column  effluent 
exceeding  95%  within  five  pore  volumes  (Fig.  4).  This 
is  consistent  with  the  minimal  adsorption  of  As04  onto 
quartz  observed  in  batch  studies  at  pH  values  >3  (Xu 
et  al.,  1991).  The  calculated  sorbed  As04  site  densities 
per  mole  free  Fe  (from  As04  sorption  maxima)  were 
similar  to  site  densities  obtained  for  Fe  oxides  in  other 
studies,  suggesting  that  the  free  Fe  oxides  present  in 
the  sand  were  responsible  for  the  retardation  and  slow 
desorption  of  As04  during  transport  at  pH  4.5  and  6.5. 

The  enhanced  transport  of  As04  at  pH  8.5  is  consis- 
tent with  the  pH  dependence  of  surface  complexation 
reactions  describing  the  adsorption  of  As04  by  metal 
oxide  surfaces,  where  sorption  of  As04  decreases  dra- 
matically at  pH  values  >8  (Xu  et  al.,  1991;  Anderson 
et  al.,  1976;  Pierce  and  Moore,  1982).  Arsenate  adsorp- 
tion by  Fe  hydroxides  can  be  categorized  as  a  ligand 
exchange  mechanism,  and  it  has  been  shown  by  Ex- 
tended X-ray  Absorption  Fine  Structure  Spectroscopy 
(EXAFS)  that  As04  forms  a  combination  of  mono-  and 
bidentate  complexes  with  surface  Fe  sites  (Waychunas 
et  al.,  1993).  For  AsQ4  adsorption  onto  goethite,  these 
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lujand  exchange  reactions  have  been  represented  by 
(Goldberg,  1986;  Grossl  and  Sparks,  1995): 

FeOH(s)  +  H,As04(aq)  =  FeH2AsO„(s)  +  H20 

FeOH(s)  +  H,As04(aq)  =  FeHAs04"(s) 

+  H20  +  H  +  (aq) 

FeOH(s)  +  H,AsQ4(aq)  =  FeAsOi"(s) 


2FeOH(s) 
2FeOH(s) 


H,As04(aq) 
H3As04(aq) 


+  H20  +  2H  +  (aq) 
=  Fe2HAs04(s)  +  2H20 

=  Fe2As04"(s)  +  2H20 
+  H  +  (aq) 


The  pH  dependence  of  As04  sorption  predicted  from 
these  reactions  as  described  by  Goldberg  (1986)  and 
Grossl  and  Sparks  (1995)  supports  the  enhanced  trans- 
port at  pH  8.5  observed  in  this  study. 

The  sorption  and  subsequent  transport  behavior  of 
As04  as  a  function  of  pH  is  dependent  both  on  the 
surface  charge  of  the  sorptive  phase  and  on  the  aqueous 
speciation  of  As04.  Zero  points  of  charge  (pHZPC)  for 
amorphous  Fe(OH),  (ferrihydrite)  and  a-FeOOH  (goe- 
thite)  range  from  8.5  to  7.8  (Stumm  and  Morgan,  1981). 
It  is  expected  that  the  pHZPC  for  the  free  Fe  oxide  present 
in  the  sand  would  fall  in  this  range;  consequently,  we 
would  expect  that  Fe-oxide  surfaces  were  predomi- 
nantly positively  charged  at  pH  values  of  4.5  and  6.5 
and  predominantly  negatively  charged  at  a  pH  of  8.5. 
The  predominant  As04  species  at  pH  values  of  4.5  and 
6.5  is  H2As04~,  and  at  a  pH  of  8.5,  the  predominant 
species  is  HAsOs  .  Consequently,  at  a  pH  of  8.5  the 
combination  of  a  predominantly  negative  surface  and  a 
more  negative  As04  solution  species  resulted  in  a  dra- 
matic increase  in  the  mobility  of  As04  (Fig.  3). 

Effects  of  Carbonate  Species 

To  test  the  effects  of  HCG%/COC  ions  (present  in 
the  C02(g)  buffered  systems)  on  the  transport  of  As04 
at  pH  8.5,  column  experiments  were  performed  in  both 
C02(g)  buffered  and  N2(g)  atmospheres  (Fig.  5).  Efflu- 
ent samples  from  both  systems  were  periodically  ana- 
lyzed for  total  alkalinity  yielding  average  values  for  car- 
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velocity  =   1  cm  h" ';  background  solution  =  0.01  M  KG;  As04 
pulse  =  1  pore  volume  at  133  u.W). 


bonate  concentrations  of  1.2  mM  and  <0.1  mM  in 
C02(g)  and  N2(g)  systems,  respectively.  A  value  of  1.2 
mM  HCOi"  is  consistent  with  equilibrium  ofpCOi(g)  = 
0.0003  aim  at  pH  8.5. 

After  10  pore  volumes  of  elution,  the  percent  of  ap- 
plied As04  recovered  in  the  column  effluent  was  100% 
in  the  presence  of  C02(g),  and  >95%  in  the  N2(g)  sys- 
tem. While  a  majority  (~90%)  of  the  applied  As04 
was  eluted  within  three  pore  volumes,  tailing  was  still 
observed  in  both  columns  (Fig.  5).  The  transport  of 
As04  in  the  C02(g)  buffered  system  exhibited  a  slight 
increase  in  recovery  after  three  pore  volumes,  indicating 
that  some  competition  of  HCGT/COi"  for  As04  adsorp- 
tion sites  may  have  occurred  during  transport.  However, 
results  from  this  comparison  confirm  that  enhanced 
transport  of  As04  at  pH  8.5  (Fig.  3)  was  caused  primarily 
by  increases  in  pH  and  not  due  to  competition  from 
HCOf/CO2,    ions. 

Effects  of  Phosphate  Competition 
on  Arsenic  Transport 

At  pH  4.5,  increasing  additions  of  P04  to  column 
influent  resulted  in  increased  recoveries  of  applied  As04 
and  caused  the  As04  BTCs  to  shift  to  the  left,  indicating 
a  decrease  in  As04  retardation  (Fig.  6).  These  results 
are  consistent  with  other  studies  showing  enhanced  mo- 
bility of  As04  in  the  presence  of  P04  (Peryea,  1991). 
After  10  pore  volumes  of  elution,  the  percent  of  applied 
As04  recovered  in  the  column  effluent  was  35.6,  36.7, 
63.5,  and  88.4  for  P04  concentrations  of  0,  13.4,  134, 
and  1340  u.M,  respectively.  With  the  addition  of  1340 
u.M  P04,  the  As04  BTC  occurred  slightly  after  one  pore 
volume,  and  had  a  peak  concentration  (CIC„)  near  1.2. 

Breakthrough  curves  of  As04  in  the  presence  of  vary- 
ing P04  concentrations  at  pH  6.5  (Fig.  7)  exhibited  the 
same  general  behavior  as  As04  BTCs  at  pH  4.5,  but  in 
all  cases  the  recovery  of  applied  As04  was  greater  at 
pH  6.5  than  at  pH  4.5,  consistent  with  the  effects  of  pH 
noted  in  Fig.  3.  After  10  pore  volumes  of  elution,  the 
percent  of  applied  As04  recovered  in  the  column  efflu- 
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Fig.  6.  Transport  of  arsenate  (As04)  through  sand  at  pH  4.5  in  the 
presence  of  0,  13.4,  134,  and  1340  u.1/  PO,  (Column  conditions: 
Pore  water  velocity  =   1  cm  h   ';  background  solution  =  0.01  M 
KCI;  As04  pulse  =  1  pore  volume  at  133  \iM). 
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Fig.  7.  Transport  of  arsenate  (As<34)  through  sand  at  pH  6.5  (CO; 
bufTered)  in  the  presence  of  0, 13.4, 134,  and  1340  \lM  P04  (Column 
conditions:  Pore  water  velocity  =  1  cm  h    ';  background  solution  = 
0.01  M  KCI;  As04  pulse  =  1  pore  volume  at  133  pM ). 

ent  was  39.2,  45.0,  67.9,  and  91.6  for  the  columns  at  P04 
concentrations  of  0,  13.4,  134,  and  1340  u.A/.  respec- 
tively. Again,  with  the  addition  of  1340  \xM  P04,  the 
As04  BTC  occurred  slightly  after  one  pore  volume,  and 
had  a  C/C0  >1. 

The  observation  that  relative  C/C0  were  >1  for  As04 
BTCs  at  pH  4.5  and  6.5  in  the  presence  of  1340  u.Af 
P04,  indicates  that  As04  was  concentrated  at  the  leading 
edge  of  the  P04  BTC.  Measurement  of  P04  concentra- 
tion in  the  column  effluent  by  IC  verified  that  P04 
breakthrough  occurred  immediately  after  As04.  While 
recovery  of  applied  As04  was  dramatically  enhanced 
by  the  addition  of  1340  \j.M  P04  to  the  applied  pulse, 
the  recovery  of  applied  As04  still  did  not  exceed  92%. 
This  may  indicate  that  certain  sorption  sites  on  the  sand 
were  specific  for  As04,  or  that  desorption  of  As04  re- 
mained kinetically  limited  even  in  the  presence  of  high 
concentrations  of  P04. 

The  total  adsorption  capacity  of  a  column  at  pH  4.5 
was  estimated  from  the  batch  isotherm  data  as  being 
13.0  fjimol  (As  or  P)  per  packed  column  of  sand.  In  a 
one  pore  volume  pulse,  the  total  input  of  As04  was  2.7 
u.mol  As,  and  the  total  input  of  P04  at  molar  P/As  ratios 
of  0,  1:10,  1:1,  and  10:1  was  0,  0.27,  2.7,  and  27  u.mol  P, 
respectively.  Consequently,  the  resulting  total  column 
loading  (As04  plus  P04)  for  column  experiments  was 
2.7.  3.0,  5.4,  and  29.7  u.mol  per  pulse  at  molar  P/As 
ratios  of  0,  1:10,  1:1,  and  10:1,  respectively.  At  a  10:1 
molar  P/As  ratio,  the  total  column  adsorption  capacity 
was  exceeded  by  approximately  100%,  indicating  that 
there  should  have  been  sufficient  P04  present  in  the 
column  to  compete  for  all  adsorption  sites. 

An  additional  experiment  was  performed  at  pH  4.5 
to  evaluate  the  effects  of  P04  competition  on  As04 
transport  after  the  applied  As04  had  time  (20  pore 
volumes  =  7  d)  to  react  with  the  solid  phase  (Fig.  8). 
The  recovery  of  As04  from  a  one  pore  volume  pulse 
of  143  u.M  As04  solution  was  only  20%  following  20 
pore  volumes  of  0.01  M  KCI.  At  20  pore  volumes,  a  1420 
[xM  PQ4  solution  (10:1  P/As  molar  ratio  with  respect 
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Fig.  8.  Transport  of  arsenate  (AsO<)  through  sand  at  pll  4.5  followed 
by  continuous  addition  of  1420  \iM  P04  at  20  pore  volumes  (Col- 
umn conditions:  Pore  water  velocity    =    1   cm  h    ';  background 
solution  =  0.01  M  KCI;  As(J4  pulse  =  1  pore  volume  at  143  (xA/). 

to  the  As04  pulse)  was  applied  as  a  continuous  pulse. 
Arsenate  recovery  increased  dramatically  by  more  than 
35%  within  15  additional  pore  volumes  of  elution.  how- 
ever significant  As04  tailing  was  observed  after  the  first 
five  pore  volumes  of  P04  addition,  indicating  that  a 
large  fraction  of  the  adsorbed  As04  was  not  readily 
desorbed,  even  in  the  presence  of  excess  P04.  Further- 
more, the  percent  recovery  of  applied  As04  in  the  efflu- 
ent was  <60%  after  15  pore  volumes  of  P04  elution 
and  35  pore  volumes  total  elution.  In  contrast,  a  simulta- 
neously applied  P04  pulse  resulted  in  a  recovery  of 
>85%  within  10  pore  volumes  (Fig.  6).  This  suggests  a 
decrease  in  As04  desorption  rates  in  'aged'  columns 
compared  to  columns  receiving  simultaneous  inputs  of 
P04  and  AsO„. 

It  is  unclear  whether  the  difference  in  As04  recovery 
(in  the  presence  of  P04)  between  'aged'  and  'unaged' 
systems  was  due  to  changes  in  the  apparent  desorption 
rates  of  surface  Fe-As  complexes,  the  formation  of  As 
precipitates,  or  the  formation  of  occluded  Fe-As  surface 
complexes.  Although  no  studies  were  performed  to 
identify  possible  As04  solid  phase  formation,  calcula- 
tions of  ion  activity  products  (IAPs)  suggest  that  our 
solutions  were  well  undersaturated  with  respect  to  As04 
solid  phases.  Furthermore,  Waychunas  et  al.  (1993)  con- 
vincingly showed  that  scorodite  (FeAs04-2H:0)  forma- 
tion on  ferrihydrite  was  not  important  even  at  As04 
concentrations  of  5  mM.  Consequently,  it  is  likely  that 
the  observed  decrease  in  As04  recovery  in  the  'aged' 
system  was  due  to  a  decrease  in  the  apparent  desorption 
rates  of  surface-bound  As.  In  'aged'  systems,  sufficient 
time  may  exist  for  As04  to  diffuse  further  into  the  sorb- 
ing  matrix,  resulting  in  apparent  decreases  in  As  desorp- 
tion rates,  even  in  the  presence  of  high  P04  concen- 
trations. 

In  a  separate  study,  Darland  and  Inskeep  (1997)  dem- 
onstrated that  the  transport  of  As04  (pH  =  4.5)  through 
this  sand  exhibited  chemical  nonequilibrium  at  pore 
water  velocities  ranging  from  0.2  to  90  cm  h"1.  These 
velocities  exceeded  time  scales  necessary  for  sorption- 
desorption  equilibrium.  The  tailing  in  As04  BTCs  ob- 
served in  this  study  in  the  presence  and  absence  of  P04 
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also  reflect  the  effects  of  slow  adsorption-desorption 
kinetics  relative  to  pore  water  velocity.  For  example,  at 
slower  pore  water  velocities  (<0.2  cm  fT1),  As  recovery 
from  the  aged  column  may  have  more  closely  approxi- 
mated equilibrium  conditions,  resulting  in  a  greater  re- 
covery of  applied  As.  Consequently,  the  effects  of  As- 
P  sorption  competition  on  As04  transport  observed  in 
this  study  cannot  necessarily  be  modeled  using  sorption 
equilibria  for  AsG"4  and  P04,  but  would  require  kinetic 
approaches  relying  on  formulation  of  empirical  (or  fun- 
damental) kinetic  expressions  describing  As04  and  P04 
sorption-desorption  (Darland  and  Inskeep,  1997). 


CONCLUSIONS 

Batch  adsorption  experiments  and  column  transport 
experiments  indicated  that  As04  was  strongly  sorbed 
to  a  sand  containing  minor  amounts  of  free  Fe  oxides 
(0.01%  Fe,  w/w).  Given  that  the  primary  constituents 
present  in  the  sand  (feldspar  and  quartz)  are  not  consid- 
ered strong  sorptive  phases  for  As04,  the  free  Fe  oxides 
present  (as  determined  using  citrate-dithionite  selective 
dissolution)  were  likely  the  principal  reactive  phases  for 
As04  sorption.  Furthermore,  the  maximum  adsorption 
density  of  0.095  mol  As04  per  mol  free  Fe  determined 
from  batch  adsorption  isotherms  is  consistent  with  val- 
ues reported  for  pure  Fe-oxide  systems.  At  a  pH  above 
the  expected  pHZPC  of  the  Fe-oxyhydroxide  phase  (i.e.. 
pH  8.0),  As04  transport  through  the  aquifer  sand  in- 
creased dramatically,  where  As04  recovery  exceeded 
95%  within  10  pore  volumes.  The  effects  of  pH  on 
the  transport  of  As04  appear  consistent  with  the  pH 
dependence  of  surface  complexation  reactions  describ- 
ing AsOj  sorption  by  oxide  minerals. 

Phosphate  was  shown  to  effectively  compete  with 
As04  for  adsorption  sites  on  the  sand  in  batch  isotherms 
as  well  as  in  saturated  transport  studies.  The  P04  compe- 
tition, however,  was  not  sufficient  to  desorb  all  of  the 
applied  As04  either  in  simultaneously  applied  pulses, 
or  in  a  column  where  AsOj  was  applied  prior  to  a  con- 
centrated PO„  pulse.  When  the  applied  P04  exceeded 
the  calculated  column  adsorption  capacity  by  twofold, 
some  As04  remained  adsorbed  onto  the  sand.  Likewise, 
when  a  high  concentration  of  P04  (1420  \lM)  was  con- 
tinuously applied  to  a  column  that  had  been  previously 
spiked  with  an  As04  pulse,  recovery  of  As04  in  the 
column  effluent  did  not  exceed  60%,  even  though  the 
total  PO.,  loading  was  greater  than  the  calculated  column 
capacity  by  more  than  two  orders  of  magnitude.  These 
results  suggest  that  As04  desorption  kinetics  may  play 
an  important  role  in  the  transport  of  As04  through 
porous  media. 
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Arsenic  Solubility  and  Attenuation  in  Soils  of  the  Madison  River  Basin,  Montana: 

Impacts  of  Long-Term  Irrigation 

C.  A.  Jones,  W.  P.  Inskeep,*  J.  W.  Bauder,  and  K.  E.  Keith 


ABSTRACT 

Arsenic  concentrations  in  850  km  of  the  Madison  and  Upper  Mis- 
souri Rivers  exceed  the  Montana  State  human  health  standard  (0.24 
|i  \1 ).  In  addition,  ground  water  As  concentrations  in  the  northern 
portion  of  the  Lower  Madison  River  Valley  are  above  the  federal 
drinking  water  standard  (0.67  u.M)  and  correlate  with  high  soluble 
As  levels  in  overlying  soils.  The  objectives  of  this  study  were  to 
determine  processes  affecting  As  solubibty  in  soils  of  the  Madison 
and  Upper  Missouri  River  Basins  and  assess  potential  impacts  of 
long-term  irrigation  with  As-rich  river  water  on  As  mobility.  Sixteen 
irrigated  (211  yr  of  irrigation)  and  nonirrigated  (never  irrigated)  soils 
were  sampled  in  four  major  regions  adjacent  to  the  Madison  and 
Upper  Missouri  Rivers.  There  were  no  significant  differences  (a  = 
0.05)  in  total  or  soluble  As  levels  between  irrigated  and  nonirrigated 
soils  within  any  of  the  four  regions.  In  addition,  sorption  coefficients 
(Kt  values)  in  six  paired  soils  selected  for  additional  chemical  charac- 
terization were  not  significantly  different  (a  =  0.05)  between  irrigated 
and  nonirrigated  soils.  Ammonium  oxalate  extractable  Fe  and  Mn 
concentrations  were  found  to  be  positively  correlated  with  Kt,  and 
Olsen  P  concentration  was  found  to  be  negatively  correlated  with  h  „ 
based  on  a  multiple  linear  regression  (r!  =  0.92).  Sequential  extrac- 
tions performed  for  a  subset  of  soils  showed  that  labile  As  concentra- 
tions were  similar  between  irrigated  and  nonirrigated  soils.  Soils  that 
had  been  irrigated  for  >100  vt  had  additional  As  sorption  capacities 
of  at  least  3  mg  kg-1  in  saturated  column  studies  following  30  yT  of 
simulated  irrigation.  Based  on  this  suite  of  data,  it  was  determined 
that  irrigation  history  has  not  significantly  affected  As  solubility  or 
attenuation  capacities  in  soils  of  the  study  area. 


Geothermal  discharge  in  Yellowstone  National 
Park  provides  a  continual  source  of  As  to  the  Mad- 
ison and  Missouri  River  system  (Stauffer  and  Thomp- 
son, 1984;  Nimick  et  al.,  1998).  Concentrations  of  As 
in  850  km  of  the  Madison-Missouri  mainstem  (Fig.  1) 
exceed  the  Montana  human  health  standard  of  0.24  u.M, 
with  concentrations  approaching  5  u.Af  in  the  Madison 
River  headwaters  (Nimick  et  al.,  1998).  In  addition,  As 
levels  exceed  the  USEPA  drinking  water  standard  of 
0.67  \x.M  in  approximately  70  km2  of  ground  water  in 
the  Lower  Madison  River  Valley  (Nimick,  1998)  and 
As  in  saturated  paste  extracts  from  soils  overlying  this 
area  range  up  to  approximately  130  \iM  (Keith,  1995). 
A  previous  study  suggested  that  evapoconcentration  of 
applied  irrigation  water  was  increasing  As  concentra- 
tions in  soils,  ground  water,  and  irrigation  return  flows 
of  the  Lower  Madison  River  Valley  (Sonderegger  and 
Ohguchi,  1988).  Based  on  this  theory,  Montana  imple- 
mented controls  in  1992  to  prevent  new  irrigation  proj- 
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ects  that  would  raise  As  concentrations  in  downstream 
waters.  To  predict  the  effects  of  current  and  future  irri- 
gation projects  in  the  Madison-Missouri  River  Basin, 
it  is  necessary  to  understand  processes  controlling  As 
transport  in  both  irrigated  and  nonirrigated  soils. 

Solubility  of  As  in  soils  can  be  controlled  by  sorption- 
desorption,  precipitation-dissolution,  and  oxidation- 
reduction  reactions.  Sorption  of  arsenate  (As(V)),  the 
thermodynamically  favored  species  in  oxidized  environ- 
ments, has  been  found  to  decrease  with  increased  pH 
(Goldberg  and  Glaubig,  1988;  Darland  and  Inskeep, 
1997a;  Jones  et  al.,  1997;  Raven  et  al.,  1998).  Surface 
complexation  of  As(V)  at  Fe  oxide  surfaces  occurs  pri- 
marily as  a  bidentate  innersphere  surface  complex,  al- 
though monodentate  surface  complexes  are  also  possi- 
ble based  on  x-ray  absorption  One  structure  spec- 
troscopy (Waychunas  et  al.,  1995).  Phosphate  additions 
have  been  shown  to  substantially  mobilize  As(V)  (Dav- 
enport and  Peryea,  1991;  Peryea,  1991;  Darland  and 
Inskeep,  1997a)  due  to  competition  for  binding  sites  at 
the  oxide  surface.  At  low  redox  potentials  (<20  mV  at 
pH  7),  arsenite  (As(III))  is  the  thermodynamically  sta- 
ble species.  Many  studies  have  demonstrated  that 
As(III)  does  not  sorb  as  strongly  as  As(V)  (Frost  and 
Griffin,  1977;  Pierce  and  Moore,  1982;  Xu  et  al.,  1991), 
and  total  soluble  As  concentrations  often  increase  with 
decreased  redox  potential  (Masscheleyn  et  al.,  1991a; 
McGeehan  and  Naylor,  1994).  Precipitation-dissolution 
reactions  involving  As  are  often  less  important  than 
sorption-desorption  reactions  in  surface  soils  because 
As(V)  solid  phases  are  generally  too  soluble  to  control 
equilibrium  As(V)  concentrations  (Robins,  1981;  Sadiq 
et  al.,  1983;  Dove  and  Rimstidt,  1985;  Nishimura  et  al., 
1987).  Conversely,  As(III)-suIfide  solid  phases,  such  as 
amorphous  As2S3  and  orpiment,  are  stable  in  highly 
reduced  systems  (Webster,  1990;  Eary,  1992)  and  may 
control  As  solubility  in  saturated  soils  and  sediments. 

Irrigation  with  Madison  River  water,  containing  ap- 
proximately 0.6  to  1.3  u.M  of  As  (Nimick,  1998)  has 
the  potential  to  alter  soil  solution  As  concentrations 
through  a  variety  of  mechanisms.  For  example,  sorbed 
As  concentrations  may  increase  due  to  irrigation, 
thereby  supporting  higher  equilibrium  solution  As  con- 
centrations. Increased  pore  water  velocities  from  irriga- 
tion also  may  affect  solution  As  concentrations  due  to 


Abbreviations:  AAS,  atomic  absorption  spectrophotometry; 
HGAAS,  hydride  generation  atomic  absorption  spectrophotometry; 
AOE,  ammonium  oxalate  extractable;  I,  irrigated;  NI,  nonirrigated; 
ICP-AES,  inductively  coupled  plasma-atomic  emission  spectroscopy; 
RSD,  relative  standard  deviation;  IAP,  ion  activity  product;  SEM, 
scanning  electron  microscope;  EDAX,  energy  dispersive  analysis  of 
x-rays;  LMadN,  lower  Madison  north;  UMad,  upper  Madison;  LMadS, 
lower  Madison  south;  UMis,  Upper  Missouri;  PWV,  pore  water  ve- 
locity. 
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kinetic  limitations  of  As  sorption.  Darland  and  Inskeep 
(1997b)  found  that  nonequilibrium  conditions  existed 
at  pore  water  velocities  as  low  as  0.2  cm  h -1  in  As 
transport  studies;  therefore,  sorption  equilibrium  may 
not  be  achieved  during  irrigation  events.  Finally,  flood 
irrigation  may  decrease  redox  potentials,  favoring  the 
formation  of  As(III),  and  possibly  decreasing  As  sorp- 
tion. In  summary,  irrigation  could  have  a  variety  of 
effects  on  As  solubility  depending  on  specific  site  char- 
acteristics. Consequently,  the  objectives  of  the  present 
study  were  to:  (i)  determine  effects  of  long-term  irriga- 
tion on  accumulation  of  soil  profile  As,  (ii)  identify  soil 
characteristics  that  control  As  sorption,  (iii)  compare 
As  lability  between  irrigated  and  nonirrigated  soils,  and 
(iv)  determine  the  potential  for  continued  As  attenua- 
tion in  soils  of  the  Madison  and  Upper  Missouri 
River  Basins. 

MATERIALS  AND  METHODS 

Sample  Collection 

Soil  samples  were  collected  during  1993  and  1994  at  16 
paired  (irrigated  and  nonirrigated)  sites  adjacent  to  the  Madi- 
son and  Missouri  Rivers,  Montana  (Keith,  1995).  Sampling 
sites  were  grouped  into  four  major  regions  (Fig.  1):  Upper 
Madison  (Sites  A1-A4),  Lower  Madison  South  (Sites  B1-B4), 
Lower  Madison  North  (Sites  C1-C4),  and  Upper  Missouri 
(Sites  D1-D4).  Sampling  was  concentrated  in  the  Lower  Madi- 
son because  ground  water  in  this  vicinity  has  elevated  As 
concentrations  (0.25-2.3  \i.M;  Nimick,  1998).  Irrigated  (I)  soils 
had  been  irrigated  with  river  water  for  >20  yr,  and  at  least 
two  soils  (Bl-I  and  C2-I)  had  been  irrigated  for  approximately 
110  yr.  Nonirrigated  (NI)  soils  had  no  history  of  irrigation  and 
were  selected  based  on  textural  and  topographical  similarities 
with  the  paired  irrigated  soils.  The  paired  soils  were  from  the 
same  soil  series  and  were  <60  m  apart  at  each  sampling  site. 
Soils  were  collected  at  discrete  depth  increments  (0-5,  5-10, 
10-20  cm),  and  samples  were  stored  on  ice  during  transport 
to  the  laboratory. 

Soil  Characterization 

Soil  samples  were  passed  through  a  2-mm  sieve,  and  a 
subsample  of  the  fine  earth  fraction  was  dried  at  70°C  for  1  wk. 
Particle  size  analysis  was  performed  using  the  hydrometer 
method  (Gee  and  Bauder,  1986)  and  organic  matter  (OM) 
was  analyzed  using  a  modified  Walkley-Black  procedure  (Nel- 
son and  Sommers,  1982).  Saturated  pastes  of  air-dried  soils 
(Rhoades,  1982)  were  centrifuged  at  6000  g  for  30  min  and 
filtered  (0.2  p.)  prior  to  analysis  for  soluble  constituents.  Solu- 
ble Ca  was  analyzed  using  atomic  absorption  spectrophotome- 
try (AAS),  and  soluble  As  was  analyzed  with  continuous  flow 
hydride  generation  atomic  absorption  spectrophotometry 
(HGAAS).  Samples  for  As  analysis  were  acidified  to  3  M 
HC1,  prereduced  in  1%  KI,  and  mixed  in  a  reaction  coil  with 
0.6%  NaBH4-0.5%  NaOH.  How  rates  were  7  and  1  mL  min-1 
for  sample  and  reagent,  respectively.  Generated  arsine  was 
quantified  in  an  air-acetylene  flame  (Perkin-Elmer  3100;  Per- 
kin-Elmer,  Norwalk,  CT)  at  193.7  nm.  Total  soil  As  was  ana- 
lyzed using  a  four  acid  heat  digestion  (Welsch  et  al.,  1990), 
and  As  levels  in  the  resulting  digests  were  quantified  using 
HGAAS.  Ammonium  oxalate  extractable  (AOE)  Fe,  Mn,  and 
Al  were  analyzed  to  estimate  free  oxide  fractions  (Loeppert 
and  Inskeep,  1996),  and  Olsen  P  was  extracted  using  the 
NaHCO)  method  of  Olsen  and  Sommers  (1982). 


Townsend 


0.40  vM 


Fig.  1.  Median  surface  water  arsenic  (As)  concentrations  (Nimick  et 
al.,  1998)  and  sampling  locations  in  the  Upper  Madison  (I  I  Mad; 
Sites  A1-A4),  Lower  Madison  South  (LMadS; Sites  III -114),  Lower 
Madison  North  (LMadN;  Sites  C1-C4),  and  Upper  Missouri  ( I  !M  is. 
Sites  D1-D4).  Levels  of  As  in  the  alluvial  aquifer  underlying  the 
LMadN  exceeded  the  federal  drinking  water  standard  of  0.67  \iM 
As,  yet  were  below  this  level  in  all  other  sampled  ground  water 
(Tuck  et  al.,  1997). 


Batch  Sorption  Isotherms 

Soils  from  six  paired  sites  (Al,  Bl,  B4,  C2,  C3,  Dl)  were 
selected  for  further  analysis  based  on  soluble  and  total  As 
levels.  Depth-weighted  composites  were  formed  by  mixing  the 
0-  to  5-cm,  5-  to  10-cm,  and  10-  to  20-cm  soil  samples  in  a 
1:1:2  mass  ratio.  Batch  sorption  isotherm  experiments  for 
As(V)  were  performed  in  duplicate,  using  a  2:1  solution.solid 
ratio  (i.e.,  20  mL  0.01  M  KC1:10  g  soil).  Initial  As  concentra- 
tions were  0,  1.3,  4.0,  13,  40,  133,  and  1330  u.M,  and  each  soil 
slurry  was  shaken  on  a  horizontal  shaker  at  approximately 
120  cycles  min"1  for  96  h.  Equilibration  time  was  determined 
by  conducting  a  kinetic  sorption  study  at  13  \i.M  As  using  soils 
Bl-NI  and  C2-NI.  Samples  were  centrifuged  at  10  000  g  for 
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10  min,  the  supernatant  filtered  (0.2  |i),  and  the  filtrate  ana- 
lyzed for  As  using  HGAAS  (described  previously).  Isotherm 
results  were  fit  to  Langmuir  and  Freundlich  models,  and  re- 
values were  found  to  be  >0.999  for  all  24  samples.  Linear 
regressions  also  were  performed,  and  the  resulting  distribution 
coefficients  (A^d  values),  describing  the  linear  portion  of  As 
sorption,  were  regressed  against  pH,  O.M.,  clay  content,  AOE 
metals  (Fe,  Mn,  Al),  Olsen  P,  and  total  As. 

A  subsequent  batch  sorption  experiment  was  conducted  to 
determine  the  extent  of  As(  V)  reduction  during  the  96  h  batch 
isotherm  experiments.  Initial  As(V)  concentrations  of  0  and 
13  (jlM  were  used  for  each  of  the  12  soils.  As(V)  and  As(III) 
were  determined  using  the  following  method  modified  from 
Masscheleyn  et  al.  (1991b).  Five  milliliters  of  filtered  sample 
was  added  to  0.5  mL  of  2  M  TRIS  buffer  (pH  =  6.5).  While 
bubbling  this  buffered  solution  with  N2  (g),  a  solution  of  3% 
NaBFLyO.1%  NaOH  was  added  incrementally  (0.1  mL  over 
30  s,  3  min  wait,  0.4  mL  over  30  s,  3  min  wait)  to  liberate 
arsine  from  any  As(III)  in  the  sample.  The  original  sample 
(As(T))  and  the  speciated  sample  ( As(V))  were  both  analyzed 
for  As,  and  As(III)  concentration  was  calculated  by  difference. 
Analysis  of  duplicate  standards  made  from  NaAs02  and 
Na2HAs04  •  7H20  salts  yielded  As(III)  and  As(V)  recoveries 
of  97  ±  2%  and  103  ±  1%,  respectively. 

Sequential  Extractions 

Paired  surface  soils  (sites  Al,  C2,  and  C3),  representing 
the  full  range  of  experimentally  determined  Kd  values,  were 
subjected  to  sequential  extractions  to  estimate  the  distribution 
of  soil-bound  As  based  on  a  compilation  of  procedures  (John- 
ston and  Barnard,  1979;  Tessier  et  al.,  1979;  Ganje  and  Rains, 
1982;  Olsen  and  Sommers,  1982).  The  method  consisted  of 
sequential  extractions  in  the  following  order:  1  M  NFLC1  (ex- 
changeable), 0.1  M  NaOH  (nonoccluded  in  Al  and  Fe  oxides), 
1  M  HC1  (Ca-bound  and  occluded  in  Al  and  Fe  oxides),  and 
30%  H2O2/0.8  M  NH„OAc  (organically-bound).  Although  se- 
quential extractions  yield  only  qualitative  information  on  As 
fractionation,  this  procedure  helped  assess  relative  As  lability 
among  irrigated  and  nonirrigated  soils.  Extracted  samples 
were  centrifuged  at  10  000  g  for  10  min  and  filtered  through 
0.2-(jl  nylon  filters.  The  H202  extracts  were  analyzed  for  As 
with  inductively  coupled  plasma-atomic  emission  spectros- 
copy (ICP-AES)  and  As  in  the  remaining  extracts  was  quanti- 
fied with  HGAAS.  Total  residual  As  was  determined  using 
heat  digestion  as  described  previously.  Quality  assurance- 
quality  control  consisted  of  the  following:  all  samples  were 
extracted  in  duplicate;  percentage  of  recoveries  were  calcu- 
lated by  comparing  summed  extract  concentrations  with  total 
As;  interferences  were  tested  by  adding  a  known  quantity  of 
As  to  each  C3-I  extract  (matrix  spikes)  and  to  each  extraction 
blank  (extraction  spikes);  HGAAS  accuracy  was  assessed  with 
USEPA  certified  standards  obtained  from  Spex  Certiprep 
(Metuchen,  NJ).  The  mean  relative  standard  deviation  (RSD) 
of  duplicate  analyses  that  had  As  values  above  the  detection 
limit  was  11.4%.  Total  sequential  extraction  recoveries  for  the 
six  samples  averaged  95%,  matrix  spike  recoveries  ranged 
from  98  to  125%,  extractant  spike  recoveries  ranged  from 
96  to  109%,  and  the  average  recovery  on  duplicate  USEPA 
standards  was  99.5%. 


Analysis  for  Solid  Phases 

Ion  activity  products  (IAPs)  for  scorodite  (FeAsO,  ■  2H20) 
and  Ca3(As04)2  were  calculated  by  assuming  the  activity  of 
Fe3+  was  controlled  by  an  amorphous  Fe(OH)3(s)  (log  Kip  = 
-39.3)  and  that  I  =  0.013EC,  where  I  =  ionic  strength  in  mol 


L  '  and  EC  =  electrical  conductivity  in  dS  m-1  (Lindsay, 
1979).  Ion  activities  were  calculated  using  the  aqueous  chemi- 
cal equilibrium  model,  GEOCHEM  (Sposito  and  Mattigod, 
1979)  as  modified  by  Parker  et  al.  (1987).  Stability  constants 
for  relevant  aqueous  complexes  were  compiled  from  thermo- 
dynamic equilibrium  data  (Sadiq  and  Lindsay,  1979;  Sposito 
and  Mattigod,  1979;  Nordstrom  and  May,  1990).  Grains  of 
sample  C3-NI,  the  soil  with  the  highest  soluble  and  total  As 
concentrations  of  the  twelve  chosen  for  further  analyses,  were 
affixed  to  Al  stubs,  coated  with  C,  and  viewed  with  a  JEOL- 
6100  scanning  electron  microscope  (SEM).  The  grains  were 
mapped  for  As,  Fe,  Al,  S,  and  Ca  using  energy  dispersive 
analysis  of  x-rays  (EDAX)  and  those  grains  containing  ele- 
vated levels  of  As  were  more  closely  investigated  to  determine 
presence  of  discrete  As  solid  phases. 

Saturated  Column  Experiments 

Column  transport  experiments  were  conducted  on  compos- 
ited surface  samples  from  soils  Bl-I,  Bl-NI,  C2-I,  and  C2-NI 
to  determine  additional  As  sorption  capacities.  Both  Bl-I  and 
C2-I  had  been  flood-irrigated  from  1882  to  1976  at  approxi- 
mately 60  cm  yr"1,  and  are  currently  sprinkler-irrigated  at 
approximately  30  cm  yr"'.  Soils  were  packed  into  acrylic  col- 
umns (3  cm  diam.,  5  cm  length)  with  final  bulk  densities  of 
1.15  to  1.25  g  cm"3.  The  soils  were  saturated  from  the  bottom 
with  0.01  M  KCI,  then  treated  with  influent  (0.01  M  KC1-1.3 
\lM  As(V))  at  a  pore  water  velocity  (saturated  conditions)  of 
1.7  cm  h_l  until  the  cumulative  loading  approached  900  cm 
(approximately  30  yr  of  simulated  irrigation).  Column  effluent 
was  passed  through  0.45  \l  in-line  filters  prior  to  sample  collec- 
tion and  analysis.  Effluent  solution  As  concentrations  were 
measured  with  HGAAS  and  the  increase  in  soil  As  calculated 
by  difference. 

RESULTS  AND  DISCUSSION 

Sample  Characterization 

Average  total  As  levels  in  irrigated  (I)  soils  were  1.5 
to  3  mg  kg"1  higher  than  in  nonirrigated  (NI)  soils  in 
both  the  Upper  Madison  (UMad;  Sites  A1-A4)  and 
Lower  Madison  South  (LMadS;  Sites  B1-B4);  although 
neither  difference  was  significant  at  a  =  0.05  (Fig.  2). 
For  comparison,  total  As  levels  in  irrigated  soils  were 
estimated  to  be  1.5  to  7.5  mg  kg"1  higher  than  in  nonirri- 
gated soils  assuming  20  to  100  yr  of  irrigation  at  60  cm 
yr"1,  an  As  concentration  of  1  u.M  in  irrigation  water, 
50%  As  retention  in  the  upper  20  cm,  and  a  bulk  density 
of  1.5  g  cm"3.  The  measured  differences  in  total  As 
between  irrigated  and  nonirrigated  soils  were  within 
the  range  estimated  in  both  UMad  and  LMadS  soils, 
suggesting  some  retention  of  applied  As  in  both  regions. 
A  reversed  trend  was  observed  in  the  Lower  Madison 
North  (LMadN;  Sites  C1-C4)  where  average  total  As 
levels  were  approximately  12  mg  kg"1  lower  in  irrigated 
than  in  nonirrigated  soils.  Although  the  measured  differ- 
ence in  As  between  irrigated  and  nonirrigated  soils  in 
the  LMadN  was  not  significant  (a  =  0.05),  the  lower 
As  concentrations  in  irrigated  soils  suggest  irrigation 
may  be  leaching  As  out  of  the  soil  profile.  Total  As 
concentrations  in  soils  within  the  LMadN  were  two-  to 
three-fold  higher  than  in  the  other  three  sampling  areas, 
again  demonstrating  the  uniqueness  of  the  LMadN  re- 
gion. The  elevated  As  concentrations  within  the  LMadN 
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Table  1.  Arsenic  sorption  coefficients  (Ka)  of  six  paired  Madison 
River  Valley  soils. 


UMad  LMadS  LMadN  UMis 


LMadN 


UMis 


UMad  LMadS 

Fig.  2.  Total  and  soluble  As  concentrations  in  the  top  20  cm  of  soil 
for  the  Upper  Madison  (UMad),  Lower  Madison  South  (LMadS), 
Lower  Madison  North  (LMadN),  and  Upper  Missouri  (UMis). 
Each  value  is  an  average  of  depth-weighted  As  levels  from  four 
sampling  sites.  Bars  containing  the  same  letter  are  not  significantly 
different  (a  =  0.05).  The  lack  of  significance  between  irrigated 
and  nonirngated  As  concentrations  in  the  LMadN  region  is  due 
to  high  variability  among  nonirngated  sites. 

were  clearly  not  caused  by  irrigation,  but  instead  were 
apparently  caused  by  natural  conditions.  For  example, 
elevated  concentrations  of  As  in  localized  areas  adjacent 
to  the  Madison  and  Missouri  Rivers  have  been  attrib- 
uted to  the  presence  of  As-rich  Tertiary  volcanic  sedi- 
ments, sorption  of  river  As  to  Fe  oxides  during  flooding, 
and  upward  transport  of  reduced,  As-rich  ground  water 
(Korte  and  Fernando,  1991;  Nimick,  1998).  Total  As 
concentrations  in  the  Upper  Missouri  (UMis;  Sites  Dl- 
D4)  were  substantially  lower  than  in  the  LMadN,  and 
were  more  similar  to  As  concentrations  in  the  UMad 
and  LMadS.  Irrigated  and  nonirngated  soils  in  the  UMis 
region  had  nearly  identical  total  As  concentrations,  as 
expected  due  to  relatively  low  As  levels  (-0.4  u.Af )  in 
irrigation  water. 

Soluble  As  concentrations  followed  the  same  trends 
as  total  As  levels  among  regions  and  between  irrigated 
and  nonirrigated  soils  (Fig.  2).  Most  notably,  average 
soluble  As  levels  in  nonirrigated  soils  of  the  LMadN 
were  five-fold  higher  than  in  irrigated  soils,  suggesting 
that  irrigation  with  lower  Madison  River  water  ([As]  ~ 
0.9  \y.M)  has  decreased  naturally  occurring  high  levels 
of  As  by  dilution  and  flushing.  Average  soluble  As  con- 
centrations in  the  nonirrigated  LMadN  soils  were  ap- 
proximately 30  ixM,  more  than  50-fold  higher  than  con- 
centrations in  the  other  three  nonirrigated  areas  (Fig. 
2).  In  summary,  although  no  significant  (a  =  0.05)  differ- 
ences were  found  between  irrigated  and  nonirrigated 
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t  ht  values  are  from  linear  portions  of  duplicate  isotherms. 

soils,  UMad  and  LMadS  soils  may  have  sorbed  some 
As  from  irrigation  water,  whereas  naturally  high  levels 
of  As  in  the  LMadN  have  possibly  been  leached  and 
diluted  by  irrigation  water. 

Sorption  of  Arsenic  in  Irrigated 
vs.  Nonirrigated  Soils 

Distribution  coefficients  (Kd),  describing  the  linear 
portion  of  As  sorption  isotherms,  ranged  from  3.4  L 
kg-1  in  C2-I  to  33.2  L  kg"1  in  Al-NI  (Table  1),  indicating 
substantial  differences  in  As  sorption  characteristics 
among  soils.  Although  the  mean  KA  for  irrigated  soils 
was  somewhat  lower  than  for  nonirrigated  soils  (14.7  ± 
12.0  L  kg-1  vs.  17.4  ±  8.8  L  kg-1),  the  difference  was 
not  significant  (a  =  0.05).  Therefore,  irrigation  has  not 
significantly  lowered  the  capacity  for  additional  As  sorp- 
tion in  the  tested  soils.  A  regression  analysis  between 
K6  and  each  of  eight  variables  [pH,  OM,  clay  content, 
AOE  (ammonium  oxalate  extractable)  metals  (Fe,  Al, 
Mn),  Olsen  P,  and  total  As]  was  performed  to  determine 
which  soil  characteristics  were  correlated  with  As  sorp- 
tion. A  three  variable  multiple  linear  regression  de- 
scribed the  majority  of  variation  in  K&  values  among 
samples  (Fig.  3).  Colinearity  was  assessed  between  vari- 
ables, and  none  of  the  three  variables  were  significantly 
related  (a  =  0.05).  The  equation 

K„  =  0.015(AOE  Fe)  +  0.038(AOE  Mn) 

-  0.62(Olsen  P)  +  3.3         (r2  =  0.92) 

suggests  that  levels  of  oxide  materials  and  NaHC03- 
extractable  P,  rather  than  irrigation  history,  control  As 
sorption  in  these  soils.  It  has  been  well  documented  that 
As  sorbs  strongly  to  oxides  of  Fe  (Hingston  et  al.,  1971; 
Raven  et  al.,  1998)  and  Mn  (Oscarson  et  al.,  1981),  and 
that  P  competes  for  As  sorption  sites  (Roy  et  al.,  1986; 
Peryea,  1991;  Darland  and  Inskeep,  1997a).  Although 
irrigation  could  decrease  As  sorption  by  decreasing  lev- 
els of  Fe  and  Mn  oxides,  or  increasing  levels,  there  were 
no  significant  differences  (a  =  0.05)  between  irrigated 
and  nonirrigated  means  for  any  of  these  three  variables 
(Table  2). 

Determination  of  As(V)  and  As(III)  species  follow- 
ing the  batch  isotherm  experiment  revealed  that  >50% 
of  the  As(V)  was  reduced  to  As(III)  in  only  C2-I  and 
C3-I,  the  soils  with  the  lowest  KA  values.  This  suggests 
that  Ki  also  may  be  a  function  of  As  valence  state  due 
to  generally  accepted  differences  in  sorption  between 
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Fig.  3.  Predicted  versus  measured  K,  for  12  Madison  River  soils. 
Average  measured  Kt  values  were  calculated  from  the  linear  por- 
tion of  each  duplicate  sorption  isotherm.  The  predicted  K,  is  based 
on  a  multiple  linear  regression  with  AOE  (ammonium  oxalate 
extractable)  Fe,  AOE  Mn,  and  Olsen  P  (all  with  units  of  mg  kg  ') 
as  independent  variables. 

As(V)  and  As(III)  (Pierce  and  Moore,  1982;  Xu  et  al., 
1991);  although  it  has  recently  been  demonstrated  that 
ferrihydrite  can  sorb  similar  amounts  of  As(III)  and 
As(V)  at  high  As  surface  coverages  (Raven  et  al.,  1998). 
When  data  from  the  two  reduced  samples  were  excluded 
from  the  regression  analysis,  the  coefficients  remained 
virtually  the  same,  suggesting  the  production  of  As(III) 
in  the  two  samples  did  not  grossly  affect  the  regression 
model.  In  summary,  levels  of  oxide  materials  and 
NaHCGyexchangeable  P  apparently  account  for  the 
majority  of  the  variability  in  Kd  values  among  soils. 

Partitioning  of  Arsenic  in  Selected  Soils 

Irrigation  has  not  consistently  altered  the  partitioning 
of  As  in  three  soils  from  the  Upper  and  Lower  Madison 
regions  based  on  a  comparison  of  sequential  extraction 
data  between  irrigated  and  nonirrigated  soils  (Fig.  4). 
More  notable  differences  in  As  partitioning  and  lability 
exist  between  sites.  For  example,  the  1  M  NH4CI  and 
0.1  M  NaOH  fractions  account  for  approximately  10  mg 
kg"1  in  both  C3-I  and  C3-NI,  yet  <2.5  mg  kg"1  in  the 
other  four  soils.  These  differences  in  the  two  most  labile 
fractions  are  consistent  with  the  high  soluble  As  concen- 
trations found  in  C3-I  (12  \iM)  and  C3-NI  (20  uJW) 
compared  with  soluble  As  concentrations  of  <3  u-M 

Table  2.  Average  chemical  and  physical  characteristicst  from  the 
top  20  cm  of  six  paired  Madison  River  Valley  soils. 


Parameter 

Irrigated 

Nonirrigated 

/>-value 

Significance 

K„  L  kg  't 

14.7  ±  12.0 

17.4  ±  8^ 

0.66 

NS§ 

pH 

8.0  ±  0.1 

8.0  ±  0.2 

0.51 

NS 

OM,  % 

2.8  ±  1.9 

3.1  ±  2.7 

0.79 

NS 

Clay,  % 

20.7  ±  9.1 

17.3  ±  83 

0.56 

NS 

AOE1  Fe,  mg  kg  ' 

656  ±  479 

455  ±  162 

0.35 

NS 

AOE  Mn,  mg  kg  ' 

292  ±  163 

322  ±  179 

0.77 

NS 

AOE  Al,  mg 

kg  ' 

828  ±  246 

789  ±  277 

0.80 

NS 

Olsen  P,  mg 

kg" 

153  i  9.2 

83  ±  3.0 

0.11 

NS 

Total  As,  mg 

kg" 

14.9  ±  7.9 

20.0  ±  21.1 

0.59 

Ns 

E 
o 


c 

01 

o 

c 
o 


t  The  pH,  OM,  day,  and  total  As  values  are  depth-weighted  averages 
(Keith,  1995),  and  the  remaining  results  are  for  depth-weighted  com- 
posites. 

t  Kt  values  were  calculated  from  the  lunar  portion  of  duplicate  isotherms. 

§  NS,  Not  significant  (a  =  0.05). 

t  AOE,  Ammonium  oxalate  extractable  (Loeppert  and  Inskeep.  19%). 


A1-I        A1-NI        C24        C2-NI        C3-I        C3-NI 
Sample 

Fig.  4.  Concentrations  of  As  recovered  from  each  sequential  extract- 
ant  for  three  irrigated  (I)  and  nonirrigated  (NI)  soils.  Values  are 
averages  of  results  from  extractions  performed  on  duplicate 
samples. 

found  in  the  other  four  soils.  In  addition,  residual  frac- 
tions account  for  <20%  of  the  total  As  in  C3,  yet  ap- 
proximately 50%  of  total  As  in  Al  and  C2,  also  indicat- 
ing large  differences  in  As  lability  between  sites.  Soils 
from  C2  and  C3  are  both  from  the  LMadN  region,  again 
demonstrating  the  high  degree  of  heterogenity  in  As 
distribution  within  this  region.  In  summary,  irrigation 
history  has  not  substantially  altered  As  distribution  at 
three  selected  sites. 

Ion  activity  products  (IAPs)  calculated  for  the  upper 
three  depth  increments  (0-5,  5-10,  and  10-20  cm)  of 
each  soil  were  compared  with  Ksp  values  (Wagemann, 
1977)  for  Ca3(As04)2  and  scorodite  (FeAs04  •  2H20). 
IAPs  were  lower  than  Ksp  values  by  at  least  six  and 
nine  orders  of  magnitude  for  Ca3(As04)2  and  scorodite, 
respectively,  indicating  saturated  pastes  from  these  soils 
were  undersaturated  with  respect  to  both  solid  phases. 
Given  inconsistencies  in  literature  Kv  values  (Sadiq  and 
Lindsay,  1981),  and  uncertainty  that  true  equilibrium 
was  reached  during  the  24  h  saturated  paste  extraction,  it 
cannot  be  definitively  concluded  from  IAP  calculations 
that  metal  arsenates  did  not  exist  in  these  samples.  Al- 
though insufficient  data  were  available  to  test  for  pres- 
ence of  other  metal  arsenates,  analysis  of  C3-NI  using 
SEM-EDAX  revealed  no  solid  phases  concentrated  in 
As.  The  saturation  indices  coupled  with  the  SEM- 
EDAX  observations  suggest  that  discrete  metal  arse- 
nate phases  were  not  controlling  aqueous  As  activities. 
Therefore,  sorption-desorption  reactions  rather  than 
precipitation-dissolution  reactions  probably  control  As 
levels  in  soils  of  the  study  area. 

Saturated  Column  Experiments 

The  two  soils  (Bl  and  C3)  chosen  for  saturated  col- 
umn experiments  had  been  irrigated  for  approximately 
110  yr  with  As-rich  Madison  River  water;  yet  both  soils 
had  additional  capacities  to  retain  As.  Soil  As  concen- 
trations increased  between  3  (C2-I)  and  7.5  mg  kg"1 
(C2-NI)  following  approximately  900  cm  of  cumulative 
loading  with  a  1.33  \lM  As  influent  (Fig.  5).  Soil  C2-I 
had  the  lowest  Kd  (3.4  L  kg"1)  of  the  12  soils  selected 
for  further  analysis  (Table  1);  therefore,  3  mg  kg-1  likely 
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Fig.  5.  Increase  in  soil  arsenic  as  a  function  of  water  loading  in  satu- 
rated soil  columns  (5  cm  long,  3  cm  diam.).  Influent  (1J3  u.M  As 
in  a  0.01  M  KC1  background  solution)  was  pumped  at  a  pore  water 
velocity  of  1.7  cm  h_l  in  both  irrigated  (I)  and  nonirrigated  (NI) 
soil  columns. 

represents  the  lower  boundary  of  additional  As  reten- 
tion capacities  in  these  soils.  Increases  in  soil  As  in  the 
saturated  column  experiments  appeared  to  be  positively 
related  to  96  h  batch  KA  values,  suggesting  that  As  reten- 
tion observed  in  the  column  experiments  also  was  con- 
trolled by  levels  of  oxide  materials  and  Olsen  P.  Re- 
moval efficiencies  of  As  during  the   first  five  pore 
volumes  ranged  from  48  (C2-I)  to  85%  (Bl-I);  however, 
by  the  end  of  this  experiment,  the  soils  were  removing 
<22%  of  the  influent  As.  The  total  mass  of  As  attenu- 
ated under  field  conditions  will  be  higher  because  the 
column  experiments  were  only  performed  with  5  cm 
columns.  Differences  in  pore  water  velocity  (PWV)  be- 
tween this  experiment  (PWV  =  1.7  cm  h"1)  and  field 
conditions  also  may  have  a  large  effect  on  As  transport. 
In  a  packed  sand  column  (Darland  and  Inskeep,  1997b), 
As  recovery  in  the  effluent  after  10  pore  volumes  was 
7.2,  35.6,  and  53.3%  at  PWVs  of  0.2,  1,  and  10  cm  h  '. 
consistent  with  sorption-related  nonequilibrium.  There- 
fore, PWVs  are  expected  to  affect  As  attenuation,  and 
actual  As  retention  is  expected  to  be  higher  than  mea- 
sured due  to  anticipated  slower  average  pore  water  ve- 
locities in  the  field.  In  addition,  anticipated  flow  condi- 
tions  in   the   field   (transient,  unsaturated   flow)   are 
expected  to  result  in  greater  As  removal,  compared  to 
laboratory  column  experiments  (steady-state,  saturated 
flow).  Conversely,  there  is  a  lower  likelihood  for  prefer- 
ential flow  in  column  experiments  due  to  generally 
fewer  macropores  in  repacked  soil  columns  compared 
with  field  conditions.  Despite  differences  between  labo- 
ratory and  field  conditions,  the  column  transport  studies 
on  soils  from  sites  Bl  and  C2,  combined  with  earlier 
results  (Darland  and  Inskeep,  1997b),  support  the  con- 
tention that  As  transport  in  these  soils  is  likely  con- 
trolled by  As  sorption  to  oxide  materials  and  pore  wa- 
ter velocities. 


CONCLUSIONS 

Soils  adjacent  to  the  Madison  and  Missouri  River 
Valley  were  studied  to  determine  if  long-term  irrigation 


was  affecting  soil  profile  As  levels,  sorption  coefficients, 
lability,  or  mobility.  Total  and  soluble  As  levels  in  the 
upper  soil  profile  (top  20  cm)  were  slightly  higher  in 
irrigated  soils  than  in  nonirrigated  soils  of  the  Upper 
Madison  and  Lower  Madison  South  regions,  suggesting 
some  retention  of  applied  As.  This  observed  trend  was 
reversed  in  the  Lower  Madison  North  where  average 
soluble  and  total  As  levels  were  6-  and  1.5-fold  higher, 
respectively,  in  nonirrigated  soils  than  in  irrigated  soils. 
The  differences  imply  that  irrigation  has  lowered  natu- 
rally-occurring high  levels  of  As  in  soils  from  the  Lower 
Madison  North  through  dilution  and  flushing.  Average 
sorption  coefficients  were  16%  lower  in  irrigated  than 
in  nonirrigated  soils  although  this  difference  was  not 
significant  (a  =  0.05).  Sorption  coefficients  were  posi- 
tively correlated  with  measured  levels  of  free  Fe  and 
Mn  oxides,  and  negatively  correlated  with  NaHCOr 
extractable  P.  Therefore,  dissolution  of  oxides  or  addi- 
tional inputs  of  P  would  both  be  expected  to  decrease 
As  sorption.  Sequential  extractions  showed  no  major 
increases  in  labile  As  between  irrigated  and  nonirrigated 
soils,  consistent  with  solution  As  data.  Saturated  trans- 
port studies  demonstrated  that  surface  soils  that  have 
been  irrigated  with  As-rich  water  for  >100  yr  have  addi- 
tional As  sorption  capacities  of  at  least  3  mg  kg-1,  and 
these  capacities  are  positively  correlated  with  sorption 
coefficients.  Although  As  removal  efficiencies  fell  from 
>48%  to  <22%  following  30  yr  of  simulated  irrigation 
in  5-cm  columns,  attenuation  of  As  is  expected  to  be 
higher  in  the  field  due  to  deeper  soil  profiles  and  slower 
average  pore  water  velocities.  In  summary,  soil  profile 
As  levels  and  sorption  characteristics  have  not  been 
significantly  affected  by  up  to  110  yr  of  irrigation  with 
As-rich  Madison  and  Upper  Missouri  River  water.  Due 
to  the  large  variability  in  soil  profile  As  levels  and  sorp- 
tion capacities  between  sampling  sites,  potential  impacts 
from  irrigation  projects  should  be  assessed  on  a  site- 
specific  basis. 
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